Dynamics and transport in random quantum systems governed by strong-randomness 
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We present results on the low-frequency dynamical and transport properties of random quantum 
systems whose low temperature (T), low-energy behavior is controlled by strong disorder fixed points. 
We obtain the momentum and frequency dependent dynamic structure factor in the Random Singlet 
(RS) phases of both spin-1/2 and spin-1 random antiferromagnetic chains, as well as in the Random 
Dimer (RD) and Ising Antiferromagnetic (lAF) phases of spin-1/2 random antiferromagnetic chains. 
We show that the RS phases are unusual 'spin metals' with divergent low-frequency spin conductivity 
at T = 0, and we also follow the conductivity through novel 'metal-insulator' transitions tuned by 
the strength of dimerization or Ising anisotropy in the spin-1/2 case, and by the strength of disorder 
in the spin-1 case. We work out the average spin and energy autocorrelations in the one-dimensional 
random transverse field Ising model in the vicinity of its quantum critical point. All of the above 
calculations are valid in the frequency dominated regime lj > T, and rely on previously available 
renormalization group schemes that describe these systems in terms of the properties of certain 
strong-disorder fixed point theories. In addition, we obtain some information about the behavior 
of the dynamic structure factor and dynamical conductivity in the opposite 'hydrodynamic' regime 
Lu < T for the special case of spin-1/2 chains close to the planar limit (the quantum x-y model) 
by analyzing the corresponding quantities in an equivalent model of spinless fermions with weak 
repulsive interactions and particle-hole symmetric disorder. 

PACS numbers; 75.10.,]m, 78.70.Nx, 75.50.Ee, 71.30.+h 



I. INTRODUCTION 

Disorder effects arising from quenched randomness are 
at the heart of many interesting and novel phenomena ob- 
served in condensed matter systems — examples include 
Griffiths singularities near phase transitions in disor- 
dered magnets (and the related phenomenon of local- 
moment formation in disordered electronic systemscl), 
metakinsulator transitions in disordered electronic sys- 
temsjj and two-dimensional phenomena such as weak,lo- 
calization and the quantum Hall plateau transitions .H 

In particular, the interplay between disorder and quan- 
tum interference leads to unusual dynamics and trans- 
port in these systems. Such effects are well understood 
for disordered quantum systems in which many-body cor- 
relations are not significant (such as disordered Fermi liq- 
uids). In contrast, relatively little is reliably known about 
the effects of strong disorder in the presence of strong cor- 
relations (say, due to electron-electron interactions in an 
itinerant electronic system, or due to exchange interac- 
tions in a system with localized spin degrees of freedom). 

However, there does exist one class of systems where 
theoretical tools are available to analyze this interplay be- 
tween strong disorder, correlations and quantum fluctu- 
ations; important examples include-|ppe-dimensional ran- 
dom antiferromagnetic spin chainsQu and raxidom quan- 
tum Ising models in one and two dimcnsionsljiZI. In these 
quantum systems, it is possible to systematically treat 
disorder and correlation using a strong-disorder renor- 



malization group (RG) technique that is designed to be 
accurate when the strength of the disorder, as measured 
by the widths of the distributions of the various cou- 
plings, is large. Such a strong-disorder approach works 
in these problems because these systems, when studied 
at ever larger length scales (and correspondingly lower 
energy scales), appear more and more disordered. More 
precisely, the low-energy effective theory obtained from 
the RG has the remarkable property that the widths of 
the distributions of the various couplings in the theory 
grow rapidly as the energy cutoff is lowered; this means 
that the RG procedure gives reliable results for the effec- 
tive Hamiltonian that governs the low-energy properties 
of the system. Moreover, the extremely strong disorder 
present at low energies in the effective theory actually al- 
lows one to straightforwardly calculate some thermody- 
namic properties and ground state correlators within the 
effective theory — this is, in essence, because strong dis- 
order implies that some particular terms in the effective 
Hamiltonian dominate over all others; calculations can 
then be performed by treating these terms first and in- 
cluding the effects of the other terms perturbatively. This 
approach has been used successfully in the past to obtain 
a wealth of information about the low-temperature ther- 
modjggngii and ground state correlators in such sys- 

Here, we exploit this simplicity that emerges at strong 
disorder to obtain the first analytical results on the low- 
frequency dynamics and transport in these systems at 
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low temperature T. Most of our results are obtained for 
T = 0; these are expected to be exact at zero temperature 
in the low- frequency limit, and to remain valid at non- 
zero temperatures for low frequencies uj ^ T. Moreover, 
in certain special cases, we can also access the regime 
u; <T. 

In the remainder of this section, we introduce the var- 
ious systems that are studied in this paper, and describe 
the organization of the rest of the paper. A brief sum- 
mary pf| some of our results has already appeared else- 
where. E3 

Our focus is on three model systems. The first model 
we consider is the one-dimensional random antiferromag- 
netic XXZ spin-1/2 chain with the Hamiltonian 

HXXZ = E [J^(^^^ ■ + ! + 44+1) + -^44+1] ' (1) 
3 

where Sj are spin-1/2 operators at lattice sites j sepa- 
rated by spacing a, and both J^^ and are random 
positive exchange energies drawn from some probabil- 
ity distributions. Such a Hamiltonian describes the low- 
energy (magnetic) dynamics of. insulating antiferromag- 
netic spin-1/2 chain compoundgi2lll3 with chemical disor- 
der that affects the bond strengths. We will also consider 
chains with slightly different probability distributions of 
the even and the odd bonds and study the effects of such 
enforced dimerization. The strength of the dimerization 
in the bonds is conveniently characterized by a dimen- 
sionless parameter S defined as 

S - In Jc - In Jo ^^-j 
var(ln Jc) + var(ln Jo) 

where Jc (Jo) represents even (odd) bonds, and the over- 
line and "var" denote correspondingly the average and 
variance over the distribution of bonds. Thus, we have 
(5 > (5 < 0) if even (odd) bonds are stronger on av- 
erage. For future reference, we also introduce the basic 
length scale in this system, 



" var(ln Jg) -|- var(ln Jo) 

Detailed information about the spin dynamics in such 
systems can be obtained by inelastic neutron scatter- 
ing (INS) experiments that directly probe the frequency 
and momentum dependent dynamic structure factor 
S"l^{k,uj). At T = 0, S°''^ik,uj) has the spectral rep- 
resentation 

S'^^{k,u;) = i ^(0|r ,|m)(m|sf |0)5(c. - E„.) , (4) 

m 

where = J2j ^ {|"^)} denote the exact 

eigenstates of the system with excitation energies Em 
relative to the ground state |0). The symmetry of Hxxz 
under rotations about the z axis implies that we can re- 
strict our attention to two independent components S^^ 



and 5'"' . The same symmetry also implies that the total 
Sjot — J2j conserved — it then makes sense to talk 
of the spin transport in such a system. We characterize 
the transport of in terms of the dynamical spin con- 
ductivity (t{uj). The real part ct'^lu) of cr(cij) is defined 
by the relation P{uj) — (T'(tj)|VJfp(ijj), where P{uj) is 
the power absorbed per unit volume by the system when 
magnetic field with a uniform gradient ViJ(aj) (with the 
field H always in the z direction) oscillating at frequency 
uj is applied along the length of the chain. From stan- 
dard linear response theory, we have the following Kubo 
formula for a' (uj) at T = 0: 

'^'(^) = 4 E I ("^1 E io) - ■ (5) 

LUL/ 

m j = l 

In the above, tj — iJj'is^ sj^-^ — Sj+iSj)/2 is the cur- 
rent operator on link j that transfers one unit of the 
from one site to the next. Here and everywhere in the 
following, the frequency uj is taken positive for notational 
convenience. Note that both S"^{k,uj) and cr'(w) as de- 
fined here are self- averaging in the thermodynamic limit. 

The second model we consider is the random antiferro- 
magnetic Heisenberg spin-1 chain with the Hamiltonian 

Hsi = ^ JjS-j ■ Sj+i , (6) 
j 

where Sj are spin-1 operators on lattice sites j, and the 
Jj are random positive nearest-neighbor exchanges; ran- 
domness in the system is characterized by a width W of 
the corresponding distribution of log-exchanges In(Jj). 
As in the spin-1/2 case, we can characterize spin dy- 
namics and transport in terms of the dynamic structure 
factor and the dynamical conductivity; the definitions re- 
main the same except for the obvious replacement of all 
spin-1/2 operators with their spin-1 counterparts. Ex- 
perimental idealizations of pure Heisenberg spin-1 chains 
are known,Ej and experimental studies of systems with 
randomness have also been reported in the recent litera- 
ture.Ej We caution, however, that the degree of disorder 
needed to destroy the gapped Haldanp-phase of a pure 
spin-1 chain appears to be quite strong ,E3 and that all our 
calculations are done only in this strong-disorder regime. 

The third problem that we consider is the one- 
dimensional random transverse field Ising model 

HrtfIM = - ^ JjCr|cr|+i - ^ hjffj , (7) 

j 3 

with random ferromagnetic interactions Jj and positive 
random transverse fields hj; here aj are Pauli spin ma- 
trices. The strong-disorder RG approach, and its con- 
sequences for the low-temperature thermodynamics and 
static correlators, have befn analyzed in greatest detail 
for this particular model £3 Also, there are extensive niu 
merical results available for some dynamical properties .Ell 
This model thus serves as a benchmark to test reliability 
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of our approach to the calculation of dynamical prop- 
erties in these strong-disorder systems — we will analyze 
various average autocorrelation functions in considerable 
detail and compare our results with the earlier numerical 
work. 

The paper is organized as follows: We begin in 
Sec. H with a general discussion of the various types of 
states that we encounter in these models, along with an 
overview of our most important results for the dynamics 
and transport in various regimes; the last part of this sec- 
tion is devoted to a general outline of the basic approach 
that is used to obtain these results. The following Sec- 
tions III, IV, and present careful derivations of our 
results for the zcro-tcmperature dynamical properties of 
the three model systems that we consider, with each sec- 
tion starting with a review of the basic RG approach used 
to study the corresponding system. In Sec. Ill we evalu- 
ate the dynamic structure factor and the dynamical con- 
ductivity in the various phases of the random XXZ spin- 
1/2 chain. This is followed, in Sec. [V, by an analysis of 
the spin conductivity in the strongly-random Heisenberg 
antiferromagnetic spin-1 chains, and, in Sec. by an 
analysis of the average local dynamical properties of the 
random quantum Ising model in the vicinity of its critical 
point. Section VI is devoted to a qualitative analysis of 



the dynamical and transport properties of the XXZ spin- 
1/2 chains at non-zero temperatures in the regime lo <T, 
along with some quantitative calculations in the XX spin- 
1/2 chain that are possible in this case because of the 
mapping to free fermions. We conclude, in Sec. VII , with 
a discussion of the possible experimental tests of some of 
our predictions for the one-dimensional random-exchange 
antiferromagnetic spin chains. Some technical details, as 
well as some additional developments slightly removed 
from the main thrust of the paper, are relegated to the 
appendices. 



II. OVERVIEW 

Broadly speaking, our results are for two types of 
states. First, there are ground states governed (and 
therefore best described by some suitable strong-disorder 
RG approach) by infinite-randomness fixed points — 
examples include the random singlet states of the spin- 
1/2 antiferromagnetic chains and the critical point of the 
random transverse field Ising model. Then, there are 
the so-called "Grifhths" phases in the immediate vicinity 
of these critical states; in these phases, the low-energy 
renormalized randomness is strong, but not infinite. 

In both cases the low-energy excitations are localized, 
but with a characteristic 'localization length' — i.e. the 
'size of the excitation' — that diverges as a power of Inw 
for energy a; ^ 0. [We emphasize that this is the state- 
ment about the (rare) low-energy excitations and is in- 
deed valid in the Griffiths phases, even though in this 
case all equal-time correlators at T = indicate a finite 



localization length; for details see the main body of the 
paper.] Apart from this logarithmically divergent 'local- 
ization length', we can also define, from the integrated 
density of states for excitations up to energy a 
length = n^^^'^ that is the typical spacing between 
these excitations in d dimensions [the results we report 
here are for d — 1, but similar phases do occur for d> 10] . 

For a ground state governed by an infinite-randomness 
fixed point, L^^ diverges at low energies with the same 
power of In lo as the typical size of the excitation. This 
means a strongly-divergent density of states at low en- 
ergy, which allows the system to behave as a conductor if 
there is a conserved quantity (e.g., spin or particle num- 
ber) to be transported. In a Griffiths phase, on the other 
hand, ^ uj^^/'^ , with z a nonuniversal dynamical ex- 
ponent that varies continuously within the phase. Here, 
the low-energy excitations are rare; they are typically 
spaced by distance i^j, which diverges as a power-law at 
low energy and thus is much larger than the excitation's 
typical size, which is diverging only logarithmically. In 
the RG language, the Griffiths phases are governed by 
lines of fixed points ending in the infinite-randomness 
critical fixed point; along such a line, the dynamical ex- 
ponent z varies continuously and diverges near the criti- 
cal point. 

In terms of the original microscopic model, the low- 
lying excitations in the Griffiths phases come from re- 
gions where the local quenched random variables deviate 
strongly from their global averages. These deviations are 
such that the local averages would put that region in a 
different phase. If the system is not at a phase transi- 
tion, the probability of such a rare region occuring and 
being of linear size L behaves as e~^^'" for large L, for 
some constant ci . Such a rare region typically results in a 
low-lying mode with a sharply defined (in the sense that 
C2, introduced below, is sharply defined) characteristic 
frequency proportional to e""^^^ . This gives rise to a 
power-law low-energy density of states, with the dynam- 
ical exponent z being determined by the ratio of the con- 
stants cxjci. For a disordered Griffiths phase, the rare 
regions are finite "islands" of either an ordered phase, 
or a different disordered phase. The resulting low- lying 
excitations localized on these rare regions produce a low- 
frequency conductivity a' (iS) or scaled dynamic struc- 
ture factor u!S{k,uj) vanishing as tj^/^ at low frequencies 
(apart from possible logarithmic factors attributable to 
singular low-energy behavior of the relevant matrix ele- 
ments that may, in some cases, be sensitive to the loga- 
rithmically divergent size of the relevant excitations). 

For one-dimensional systems, there are also power-law 
Griffiths effects in Ising-ordered phases. These occur be- 
cause of rare regions locally in the disordered phase. The 
low-energy excitation associated with such a region is 
a domain wall (or "kink"). To produce a single such 
low-energy domain wall requires flipping the spontaneous 
magnetization on one side of the the wall, which is tan- 
tamount to flipping a semi-infinite piece of the chain. 
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Such a flip of an infinite domain cannot occur at a fi- 
nite (non-zero) frequency. Tlie leading contribution to 
the low-frequency dynamics is then associated with two 
nearby such rare low-energy domain walls which allow 
the ordered domain between them to flip at a low but 
nonzero frequency. The result of this is that the low- 
frequency (j'{uj) and ujS{k,Lu) vanish as w^/^ at low fre- 
quency in these one-dimensional Ising-ordered Griffiths 
phases (we are again ignoring possible logarithmic fac- 
tors which can arise for precisely the same reasons as in 
the disordered phase). Note however that the Griffiths 
singularities in Ising-ordered phases in d > 1 are of a 
very different character — in these cases, the low-energy 
density of states vanishes faster than any power oi uj, as 
is discussed in Ref. 0. 

In Sections |n|-0 we will provide a detailed justification 
of these general observations by explicitly calculating the 
low- frequency dynamical properties in a variety of cases. 
In the rest of this Section, we review the phase diagrams 
of our model systems, and highlight our most important 
results in each case. 



A. Random antiferromagnetic XXZ spin-1/2 chains 

1. Phase dianram 



The phase diagram of the random antiferromagnetic 
XXZ spin-1/2 chains is best understood as a product of 
the competition between the transverse part of the cou- 
pling J-^, which favors singlet formation, and the 'clas- 
sical' interaction J^, which favors a ground state with 
Ising antiferromagnetic order. 

When the dominate, the ground state can be 
loosely thought of as being made up of singlet pairs. In 
this Random Singlet (RS) state, the interplay of disorder 
and quantum fluctuations locks each spin into a singlet 
pair with another spin; the two spins in a given singlet 
pair can have arbitrarily large spatial separation, with 
the disorder determining the particular pattern of the sin- 
glet bonds in a given sample. On the other hand, when 
the dominate, the system has Ising antiferromagnetic 
(lAF) order in the ground state (with the spins all ori- 
ented parallel to the z axis), although Grifhths effects can 
fill in the gap leading to an lAF ordered Griffiths phase. 

These two states are separated by a quantum phase 
transition that occurs when the couplings and 
have roughly similar distributions (have roughly equal 
strengths). A special feature of this system is that the 
ground state at any point on the critical manifold is also 
a Random Singlet state, though the details of the exci- 
tation spectrum are somewhat different. 

If we now turn on enforced bond dimerization starting 
with the RS state that obtains for small J^, or the RS 
state of the Heisenberg chain, the system moves into a 
Griffiths phase dubbed the Random Dimer (RD) phase; 
in this phase the singlet bonds in the ground state now 




AF: 



FIG. 1. Schematic phase diagram of random antiferromag- 
netic XXZ spin-1/2 chains obtained in Ref. ^, showing the 
three different random singlet (RS) fixed points and RG fiows. 
Our prediction for the low-frequency behavior of the dynam- 
ical conductivity is indicated for each phase. For details, see 
here z = z{Siaf) is a (continuously varying) dynam- 
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ical exponent in the lAF Griffiths phase. 
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FIG. 2. The Random Dimer phases for XX or Heisenberg 
spin-1/2 antiferromagnetic chains, represented as lines of fixed 
points ending in the critical fixed point labeled RS that de- 
scribes the Random Singlet state at zero dimerization;EI here 
z — z{5) is a dynamical exponent in the RD Griffiths phase. 



preferentially start on one sublattice and end on the 
other. 

Schematic phase diagrams summarizing the above are 
shown in Figs. |l| and |[ 



2. Spin transport 

We characterize the spin transport properties of the 
various phases in terms of the low-frequency behavior 
of the dynamical conductivity: We find that the T — 
dynamical conductivity diverges at low frequencies in the 
RS phase as well as at the RS critical points as 



a'{uj) = ICjislv^Lo, 



(8) 



where we have taken the opportunity to introduce the 
log-energy scale 



ln(r2o/^)- 



(9) 



Here and henceforth we use fig to denote the non- 
universal microscopic energy cutoff, which corresponds 
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roughly to the energy scale in the bare Hamiltonian for 
our various models; also, we use ly to denote the non- 
universal microscopic length scale in the problem. For 
the XXZ spin- 1/2 system near the RS phase and with suf- 
ficiently strong disorder, which is what we assume in the 
following, the microscopic length ly is given by Eq. (|^). 
[If, on the other hand, the bare disorder is weak and the 
system flows to strong disorder, then ly is the length scale 
at which the streiigth of the disorder becomes of order 
one.] /Crs in Eq. (g) is an order one numerical constant. 
The RS phase and the RS critical points separating it 
from the lAF phase are thus unusual spin conductors. 

On the other hand, the lAF Griffiths phase is a spin 
insulator with the low-frequency T = dynamical con- 
ductivity 

•'(c^) = ICiafIv {uj/nof^''^^ \n^{no/uj) , (10) 

where ziaf('^iaf) is a (continuously varying) dynami- 
cal exponent diverging at the critical point as ziaf ^ 
'^lAF ^"^^^i and /CiAF is a non-universal amplitude van- 
ishing at the transition as /Ciaf ^ '^iaf'^'^^- Here we 
parametrized the distance from the transition to the RS 
phase by 5iaf = A — Ac (where A = /J-^). The ex- 
ponent A is the relevant RG eigenvalue controlling the 
fiow away from the critical fixed point describing the 
generic transition between the RS phase and the lAF 
phase, and the exponent ■0 characterizes the low-energy 
spectrum above the RS gr ound state at this critical point 
(see Ref. || and Sec. [II A for details). The above result 
is expected to hold in the frequency regime uj ^ ^Siaf 
with the cross-over scale ^^Sjaf given in terms of Siaf as 



cr' 



-(2-V)/A 
lAF 



ln(OoM:.J~^, 

Similarly, the RD phases are also spin insulators^ with 
the T — low-frequency dynamical conductivity 



a'iuj) = /CRD/.(^/r!o)'/'"° ln{no/Lu) ; 



(11) 



the dynamical exponent zrd (S) in the RD phase diverges 
at the transition as zrd ~ and the non- universal 

amplitude /Crd vanishes at the transition as /Crd ^ \S\. 
As in the lAF phase, this result is valid at frequencies 
well below the corresponding crossover scale fig (which 
can be also viewed as the conductivity pseudo-gap scale) ; 
in the RD phases ln{no/^s) ~ 1/|'5|. 

Thus, in both the lAF phase and the RD phase, the 
conductivity has the functional form 



(12) 



with the non-universal exponent a vanishing at the cor- 
responding transition. Note that a similar form but with 
fixed a = 2 — the Mott formula — is obtained via the 
usual Mott argument for the T — dynamical conductiv- 
ity of the one-dimensional Anderson insulator (the fixed 
value of a in this case simply refiects the fact that low- 
energy density of states in the Anderson insulator is con- 
stant, in contrast to the situation in the Griffiths phases 
of interest to us here). 
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FIG. 3. Sketch of the dynamic structure factor at fixed 
cj <C Slo in the RS states. 
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FIG. 4. Sketch of the dynamic structure factor at fixed 
<C f2« in the RD phases. 



3. Spin dynamics 

Turning to the spin dynamics, we find that the T — 
dynamic structure factor in the RS states in the vicinity 
of fc = 7r/a can be written in the following unusual scaling 
form 



a 



A 



^(\qly\'/''lliino/Lo: 



(13) 



for |g| <C a^^ and uj ^ flo; here a/? = H — or zz, A is 
an order one numerical constant, ly is the microscopic 
length defined earlier, and the fully universal function 
<I>(x) 



calculated in Sec. Ill 



interpolates smoothly be- 
tween the limiting forms ~¥lx) ~ 2 — 7x'^/90 for x ^ 1 
and $(x) ~ 1 -I- 2a;(cosa; -|- sina;)e^^ for a; ^ 1. A plot 
of the momentum dependence of the dynamic structure 
factor near k = ir/a (at fixed low frequency) is shown in 
Fig. H; an interesting aspect is the non-monotonic nature 
of the lineshape. We will see in Sec. [II that this oscil- 



latory behavior becomes more pronounced and leads to 
really striking structure in the momentum dependence 
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of the dynamic structure factor at (fixed) low frequency 
a; <C in the Random Dimer phases — a plot of the ex- 
pected k dependence is shown in Fig. ^. A very similar 
dependence is also predicted in the lAF Griffiths phase 
close to the transition to the RS state. 

As mentioned earlier, these results are expected to re- 
main valid at small non-zero temperatures so long as the 
frequency to satisfies uj ^ T. In Sec. VI, we will see 



that we can partially overcome even this restriction in 
the vicinity of the XX point. 



B. Spin-1 Heisenberg antiferromagnetic chains 

1. Phases 

The effect of randomness on antiferromagnetic Heisen- 
berg spin-1 chains is even more interesting. For spin- 
1, there are, in general three distinct phases possible in 
the presence of disorder. If the disorder is weak, and 
the support of the probability distribution P{J) of the 
exchanges is confined to a narrow enough region near 
the mean, then the system remains in the usual gapped, 
topologically ordered Haldane state. For stronger dis- 
order, or when P(J) has tails to large or small enough 
J, one has the 'Gapless Haldane' (GH) phase in which 
the system still has the topological order that character- 
izes the Haldane state, but becomes gapless due to Grif- 
fiths effects. Finally, if the disorder is extremely strong, 
with the (bare) distribution of exchanges broad on a log- 
arithmic scale, a Random Singlet state completely anal- 
ogous to the one encountered in the spin- 1/2 case ob- 
tains. While the GH state and the RS state are sepa- 
rated by a quantum critical point with universal critical 
properties (these propertie£u,are in fact controlled by a 
strong-disorder fixed pointO'Ll) , the corresponding transi- 
tion between the gapped and gapless Haldane states is 
a non-universal feature of the phase diagram, depending 
sensitively on the nature of the initial distribution of cou- 
plings (see Fig. ^ for a summary of the universal aspects 
of the phase diagram). 



Gapless Haldane 

XX XX XX XX 



c Random Singlet 
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FIG. 5. Schematic phase diagram of the strongly disor- 
dered Heisenberg spin-1 chain ,□ along with our results for 
the dynamical conductivity in various regimes. Moving to 
the right along the horizontal axis corresponds to increasing 
disorder. 

However, it is still possible to calculate transport prop- 
erties, such as the dynamical conductivity, that are in- 
sensitive to the details of the spatial structure (this is, 
in essence, a consequence of spin conservation). At the 
critical point separating the Gapless Haldane state from 
the Random Singlet state, we find for lo <ti ^Iq 



(j'{uj) = /Chy'd \ii'{VlQ/uo) , 



(14) 



which is a stronger divergence than in the strong-disorder 
RS phase; here ly is the non-universal microscopic length 
scale beyond which the effective model applies, ilo is the 
corresponding microscopic energy scale, and /Cry is an 
order one numerical constant. Thus, the critical point 
separating the RS phase from the GH phase is also an 
unconventional 'spin metal'. The GH phase, on the other 
hand, is a 'spin insulator', not unlike the RD phase of 
spin-1/2 chains. We find for the conductivity in the GH 
phase 



a'{uj) ^ JCGnU^/^of '''''' \n\no/u) 



(15) 



The dynamical exponent zqh varies continuously in the 
Gapless Haldane phase, diverging at the critical point 
as zgh ~ {Wc — W)^'^^^, while the the non-universal 
amplitude /Cgh(T^) remains non-zero as one approaches 
the critical point. In the above, Wc is the critical value 
of the bare disorder (the parameter W has already been 
defined in Sec. |), and the correlation length exponent 
V = 6/ (yTS — 1) is known from the analysis of Refs. |^ 
and \m. 



2. Overview of results 

In the spin-1 RS state, we obtain the same results for 
the dynamic structure factor and spin conductivity as in 
the spin-1/2 RS state, as the low-energy behavior of the 
RS state does not depend on the spin magnitude except 
through the values of some microscopic scale factors. Un- 
fortunately, once we move away from the Random Sin- 
glet state, it is difficult to discuss reliably the momen- 
tum dependence of the dynamic structure factor of the 
original spin-1 chain — this is because our actu al calc ula- 
tionsj^e done in an effective model (see section IV A and 
RcfsE'El for details) in which much of the spatial infor- 
mation about the original system is missing. 



C. Random quantum Ising spin chains 

1. Phases 

The self-dual nature of the random transverse field 
Ising model in one dimension implies that the system 
will be in a critical state if the distributions of bonds and 
fields are identical. The deviation from criticality may 
be parametrized by 



5 = 



In /i — In J 



var(ln h) + var(ln J) ' 



(16) 



with 5 > corresponding to the quantum disordered 
paramagnet, and 5 < Q corresponding to the ordered fer- 
romagnet. (Note that we use '5' both as a dimensionless 
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measure of dimerization in spin-1/2 chains, and in the 
present context; there is however no cause for confusion 
and the meaning will always be clear from the context in 
what follows.) 

This quantum critical point is flanked, for small \S\ on 
either side, by paramagnetic and ferromagnetic Griffiths 
phases with gapless excitations. 



2. Overview of results 

As mentioned earlier, these Griffiths phases and the 
quantum critical point separating them are among the 
best understood examples of such strong-randomness 
phenomena. However, all previous analyses of the dy- 
namical properties relied on numerical results supple- 
mented by scaling ideas. 

In contrast, our approach allows us to analytically cal- 
culate the average local autocorrelations of both the spin 
and the energy operators at, and in the vicinity of, the 
quantum critical point, as well as obtain the scaling be- 
havior of the dynamic structure factor of the spins. The 
main features of the average autocorrelations (as well as 
distributions of autocorrelations, which we do not ad- 
dress here) have already been noted in the earlier nu- 
merical work (Refs. |l^), while our results on the dynam- 
ical structure factor are new. Here, we only highlight 
some of the subtleties, missed in these numerical stud- 
ies, that our analytical work has uncovered regarding the 
autocorrelations — a complete tabulation of our predic- 
tions (and their interpretation in terms of Griffiths ef- 
fects) is given in Sec. 0. 

Our results for the T = imaginary-time off-critical 
spin autocorrelation in the bulk have the form 



[Clod, ir) 



|lnr| 

T"/^(<5) ' 



(17) 



where z{d) is the continuously varying dynamical expo- 
nent characterizing the Griffiths phases (from the results 
of Ref. § z-i « 2\S\, for small enough 5). In the above, 
the parameter n distinguishes between the disordered 
and ordered phases with n = 1 in the disordered phase 
and n = 2 in the ordered phase. Thus, the exponent 
controlling the power-law decay in the ordered Griffiths 
phase is twice z"^, while the corresponding exponent in 
the paramagnetic Griffiths phase is z~-^. This reflects the 
physical distinction between the disordered and the Ising 
ordered Griffiths phases noted in our general discussion 
at the beginning of this Overview. Moreover, the auto- 
correlations in the Griffiths phases are not pure power- 
law, but have a logarithmic correction, which reflects the 
fact that the appropriate 'spin' degrees of freedom rele- 
vant at a time-scale r have an effective moment of order 
I Inrj. Both these subtleties have been ignored when ex- 
tracting the dynamical exponent from the numerical re- 
sults for the average spin autocorrelations via the ansatz 



locj 



ri/^(5) ' 



and this could account for some 



of the discrepancies observed in the numerical studies. 
Similar remarks apply to other average autocorrelations 
considered, and we refer to Sec. ^ for details. 



D. The basic strategy 

We conclude with an overview of the basic strategy 
introduced by us in Ref. |l^ for the calculation of dy- 
namical and transport properties — we will be using this 
approach over and over again in what follows, and while 
the details will differ from calculation to calculation, the 
basic approach will remain unchanged. 

Consider, for concreteness, the calculation of the dy- 
namic structure factor S"^{k,uj) for the Hamiltonian 
'Hxxz- The basic idea is to eliminate high-energy de- 
grees of freedom using an appropriate strong-disorder 
renormalization group procedure (in this case, the sin- 
glet RG reviewed in Sec. HI A ), and trade in the spectral 
sum Eq. (|^) for a sum over the eigenstates of the renor- 
malized Hamiltonian Tixxz, which has fewer degrees of 
freedom and renormalized bond strengths. This renor- 
malized spectral sum must use the matrix elements of the 
renormalized versions of the spin operators; these renor- 
malized operators are of course defined by the require- 
ment that their matrix elements between the eigenstates 
of the renormalized problem reproduce the matrix ele- 
ments of the original operators between the correspond- 
ing eigenstates of the original problem. In the systems 
of interest to us, the low-energy renormalized random- 
ness is very large. In the renormalized problem at the 
energy cutoff ft ^ Hq, the effective bonds thus have a 
very broad distribution characteristic of the fixed point 
to which the system flows in the low energy limit. This 
allows us to reason as follows: Focus on pairs of spins 
coupled by 'strong' bonds in the renormalized problem, 
with strengths equal to the cutoff il. The broad distri- 
bution of bonds implies that these pairs are effectively 
isolated from their neighbors. It is therefore possible to 
unambiguously identify the excited states of these pairs 
with excitations of the full system at the same energies 
and work out the matrix elements connecting these to the 
ground state using the renormalized operators. Thus, 
to calculate the spectral sum Eq. (|), the RG IS run 
till the cutoff f2 equals rJflnai, and the problem is re- 
duced to calculating the renormalized spectral sum in 
this new theory; f2finai is chosen so that the energy of 
such excited states (associated with these strong bonds) 
measured from the ground state equals lu. The calcula- 
tion of S°'^{k,uj) then becomes a counting problem. One 
uses the known statistical properties of the renormalized 
bonds in the theory with cutoff Jlfinai to calculate the 
number of such strong bonds, and simply adds up their 
contributions weighted by the corresponding matrix el- 
ements to obtain the required result. This result is ex- 
pected to be asymptotically accurate in the limit of small 
LU, since these contributions clearly dominate in the low- 
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frequency limit. A certain simplicity thus emerges when 
the low-energy effective theory has strong disorder, and 
we will exploit this to the fullest in what follows. 

III. DYNAMICS AND TRANSPORT IN THE 
S=l/2 XXZ CHAINS 

A. Detailed characterization of the phases 

1. Singlet RG description of the Random Singlet states: A 
review 

We begin by notingi-Lhat the weak-randomness analy- 
sis of Doty and FisherEj implies that randomness is rele- 
vant for pure antiferromagnetic XXZ spin-1/2 chains for 
< J^/J"*" < 1; any amount of randomness is thus ex- 
pected to drive the system to strong disorder in this entire 
regime. 

In the stafflg-disorder regime, the singlet RG proceeds 
as follows:E2lQ We look for the bond with the largest J-^ 
coupling, say between spins 2 and 3; this sets the 
energy cutoff = maxfj^*-}. We first solve the corre- 
sponding two-spin problem and introduce the neighbor- 
ing bonds later as a perturbation. So long as the 
couplings are not large compared to the couplings, 
the ground state of the two-spin problem will always be 
a singlet separated by a large gap from the triplet ex- 
cited states. We can then trade our original Hamiltonian 
in for another Hamiltonian (determined perturbatively 
in the ratio of the neighboring bonds to the strongest 
bond) which acts on a truncated Hilbert space with the 
two sites connected by the 'strong' bond removed. To 
leading order, this procedure renormalizes the Hamilto- 
nian ?^4sitcs = Yl]=l[Jj' ^^1^1+1 + + '^/'^i^l-Hi] 
7ii4 = j^[slsl + s\sl) + Jlslsl with JjL = J^J^/{Ji + 
J|) and Jf = Jf J|/2 J;^; note that the length of this new 
bond is li = h + h + h- This procedure, if it remains 
valid upon iteration, thus ultimately leads to the Ran- 
dom Singlet state described in the overview. 

A complete understanding of the possible phases then 
requires an analysis of the effects of iterating the basic 
RG procedure. Such an analysis was performed in Ref. ^ 
leading to the following conclusions (see Fig. |l|): So long 
as the couplings do not dominate over the J-^ cou- 
plings and therefore do not produce a state with Ising an- 
tiferromagnetic order, the ground state is a Random Sin- 
glet state. In this detailed characterization of the 
low-energy effective Hamiltonian is best couched in terms 
of logarithmic variables as follows: Let fl = max{ Jj'-} 
at any given stage of the RG, and define the log-cutoff 
r = \n{ilo/^)- Also define log-couphngs Q = \n{il/jj-) 
and log-anisotropy parameters Dj = In(Aj), where Aj = 
Jj/Jj-. As r increases, the fraction of remaining sites 
np at log-cutoff scale T is given as np ^ l/T^. When 
the J"*" couplings dominate, the system rapidly flows 



to the 'XX-RS' fixed point and the probability distri- 
bution P(C, A, l\T) that determines the strengths and 
lengths of the bonds connecting the remaining sites in 
the effective Hamiltonian quickly converges to the fol- 
lowing scaling form characteristic of the XX-RS fixed 
point: P(C, A, ;|r) = Y^Vii^, j^) X ^(A). The func- 
tion Vi has been characterized in detail in Ref. ^ here 
we only note that j dyVi{x,y) = e~^. Between the lAF 
phase and this XX-RS phase lie two kinds of critical 
points. If the initial problem has full Heisenberg sym- 
metry (J^ — J-*- for each bond), the low-energy effec- 
tive Hamiltonian preserves this symmetry and has bond 
strengths and lengths drawn from the same probabil- 
ity distribution: P(C, ^ | F) = y^Vi{^, -^). In the RG 
language, the Heisenberg system is critical and is con- 
trolled by the 'XXX- RS' critical fixed point. Finally, 
in this language, the generic critical point between the 
lAF phase and the XX-RS phase is controlled by the 
'XXZC-RS' fixed point — the low-energy effective theory 
has bond strengths and lengths drawn from a distribu- 
tion P(C, DJ\V) = r3Wp2(f, ^) with V < 1 
and Jdy'P2(x,y, z) — Vi{x,z). Notice that these scal- 
ing forms imply that the distributions of the couplings 
become infinitely broad as 17 ^ 0; thus, the RG becomes 
asymptotically exact at low-energies and, in particular, 
predicts the ground state properties and low-temperature 
thermodynamics correctly. 

2. Scaling Description of the Ising Antiferromagnet 

On the Ising Antiferromagnet side, the singlet RG 
becomes invalid at low energies, and the system has a 
ground state with lAF order. The proper characteriza- 
tion of the system at these low energies is in terms of 
lAF-ordered spin clusters, as well as the domain- wall ex- 
citations that act to disrupt this order. This section is 
devoted to providing such a description. In what follows, 
we will be considering mainly the lAF phase close to the 
transition to the RS state. In this regime, the system 
will 'look' lAF ordered only well below a crossover en- 
ergy rJ^iAF I while resembling a critical system controlled 
by the XXZC critical point above the crossover scale. 
^SiAF is the scale at which the singlet RG breaks down 
and is thus determined by the properties of the RG fiows 
in the vicinity of the XXZC critical point. The corre- 
sponding log-energy scale F^jaf = Inifla/rtsiAF) is given 
ag3 TsiAp ~ ^LAF' "^itii ^ = (2 — V')/-^! where A is the lead- 
ing relevant RG eigenvalue at the XXZC fixed point and 
-0 has already been defined in the previous section. Be- 
low, we construct a scaling description of the lAF phase 
near criticality by combining information obtained from 
the singlet RG about the nature of the system at this 
crossover scale, with a 'cluster RG' approach that is de- 
signed to work in the limit of low-energies (well below 
ri^jAp) above the lAF ordered ground state. 

We begin with a sketch of our cluster RG approach. 
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Consider the Hamiltonian Tixxz with couphngs com- 
pletely dominating the couplings. Now, spins tend to 
order antiferromagnetically, and we can try formulating 
a cluster RG similar to that for the ordered phase of the 
random transverse field Ising model. Consider combin- 
ing two such spins, say S2 and S3, coupled by a strong 
bond Jl into a new 'superspin' S(23). If we identify the 
two states 1 11(23)) and | -11(23)) of this superspin with the 
states I t2i3) and 1 12^3) (which is not a unique choice), 
and treat the J-*- couplings to second order in pertur- 
bation theory, the effective Hamiltonian that we obtain 
is 

Wl(23)4 = -^1515(23) " "^3^4^(23) + '^(23)S(i23) ~ 
~ '^1(23)4(^1^ 5(23)^4 +^1^(23)^4)) 

where h,^23)=Jt, Jt(2Z)4:^'^iJt I Ji=Ji+{Jif I J2, 
and J§ = J§ + ( ^3")^/ ^2 • Thus, we see that new terms, 
not present in the original Hamiltonian, are generated: 
an effective transverse field, which acts to flip the new 
spin, and also a three-spin exchange interaction. Be- 
fore we proceed, a couple of comments regarding the new 
terms: The effective transverse field appears because the 
ground state of 7^23 is not exactly a degenerate doublet 
(the two lowest eigenstates, which are the symmetric and 
antisymmetric combinations of | T2i3) and | i2T3)) are ac- 
tually split by a small energy J^). Note also that the 
three-spin term does not violate spin conservation — for 
example, if we consider coupling the conserved total s^ot 
to a magnetic field, we immediately realize that the su- 
perspin S(23) does not couple to this field. 

In principle, we may proceed with such a clustering 
process, keeping track of all additional one- or two- or 
multi-spin- flip terms that are generated. While this RG 
is not analytically tractable, we do not expect the gener- 
ated terms to have any drastic consequences, since they 
generally become weaker and weaker, while the cou- 
plings remain almost unchanged. Alternatively, we can 
remedy this proliferation of new couplings by combin- 
ing an odd number of spins at a time — because of the 
symmetries of the Hamiltonian, any odd length chain 
will have a degenerate pair of ground states with the 
total s^ot = ±i. In addition, three-spin terms of the 
form encountered previously will now be forbidden by 
spin conservation. More explicitly, if we combine three 
spins, say S2, S3, and 54, with relatively strong couplings 
Jl and J|, into a new superspin S2 = 5(234), and treat 
the J"*" couplings perturbatively, the XXZ form of the 
effective Hamiltonian is preserved, with the new cou- 
phngs Jf2 = Jf + {J^r/{2Ji), Ji, = Ji + {Jtf/{2Ji), 
Ji2 = —2Ji J3 /J|, and J25 = — 2J2 J4 /J§. 

Either way, we will have effective spin-half objects with 
dominant Ising AF interactions. Almost always, we will 
be decimating strong couplings making larger and 
larger clusters, with the other couplings remaining 
essentially unchanged, and the remaining J-^ couplings 
growing weaker and weaker. Only rarely will there be a 



bond with the coupling large compared to the neigh- 
boring couplings, and this will then produce a singlet. 
Thus, the picture that emerges is very reminiscent of the 
ordered phase in the RTFIM. 

We may now combine this schematic cluster RG de- 
scription valid at low energies, with information about 
the crossover region obtainable from the singlet RG. At 
the crossover scale, the distribution of = \n{nsjj^p / J^) 
is given as PiC\Ts,,^p) ^ Tg^l^ exp{-C Z'^Sjaf)- Roughly 
speaking, beyond the crossover scale, the cluster RG 
merely eliminates the strongest bonds from this distri- 
bution, but keeps the low-energy tail of the distribu- 
tion unchanged. We thus expect a line of (classical) 
lAF fixed points, with properties varying smoothly with 
the distance from the criticality. The density of spin 
degrees of freedom nr in the renormalized theory is ex- 
pected to decrease as nr ~ ^s^p^~''^^^^^^^ below the 
crossover scale fi^j^p, with c some order one constant. 
This immediately gives us the density of states p{uj) ~ 
w'^nr^r^j^p ~ SflpLO~'^~^''^^''^^^ , with the continuously 

varying dynamical exponent Ziaf ^ ^iaf- The typical 
size of the excitations dominating the density of states 
scales as l^wi^) ^ Iv^Si^f^u and is much smaller than 
their typical separation ~ o;"^/^^*^. This can be readily 
seen from the qualitative picture of 'preformed tails': the 
length ldw{^) of a renormalized bond with = in the 
theory with cutoff co <C ^^^iaf scales in the same way as 
the length of a bond with ln(05,^p/a;) in the theory 
at the crossover scale f^ijAp- On the other hand, the dis- 
tribution of the log-couplings C""" = ln(ri/J-'-) is expected 
to broaden exponentially as a function of T: for example, 
when we combine n spins that are active at the crossover 
scale into a cluster, the effective transverse coupling act- 
ing on this cluster will be of order C,-^ ~ nF^j.^p. This 
then is our scaling picture for the lAF phase; the im- 
portant conclusion that emerges from this analysis is the 
fact that the transition to the RS state is preceded by a 
Griffiths phase — the lAF Griffiths phase — with a contin- 
uously varying power-law singularity in the low-energy 
density of states. 

Finally, it is also possible to obtain a rather direct 
identification of the low-energy modes in terms of the 
rare regions that dominate the low-energy dynamics; we 
conclude by sketching this briefiy here. The RG picture 
suggests that in the lAF phase the typical excitations 
at low energies w are classical domain wall excitations 
that live on the effective bonds with weak effective cou- 
plings ^ Lo. Such a weak effective can appear 
only across a long region that is locally in the RS phase. 
More quantitatively, a region of length L locally in the 
RS phase effectively corresponds to a weak bond with 
~ i7oe"'^^^. The number density of such regions is 
roughly ~ for some p <1. The density of such regions 
with < io is thus some power of u) that we choose to 
write as n{bj) ^ uj^/^- with some exponent z; the most nu- 
merous such regions will have some "optimal" (for a given 
(5iaf) microscopic structure, but whatever this structure 
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is, the corresponding optimal exponent z can be directly 
identified with the dynamical exponent 2;((5iaf) of this 
phase. This picture thus predicts that the typical separa- 
tion of such regions is of order a;~^/^, while their lengths 
are only of order | lna;|, in complete agreement with the 
schematic RG approach. 



3. Singlet RG description of the Random Dimer phases: A 
review 

While the effects of dimerization are not understood 
in detail in all regimes, it is possiblecl to use the singlet 
RG and follow the flows for a chain with full Heisenberg 
symmetry and for a chain in the vicinity of the XX-RS 
point. In these cases, a mapping to the off-critical flows 
of the RTFIM provides a detailed characterization of the 
so-called 'Random Dimer' (RD) phases that result. In 
either case_the picture that emerges can be summarized 
as follows For concreteness, assume S > 0. If disor- 
der is strong and 6 <^ 1, then the even and the odd 
bonds renormalize essentially as in the corresponding RS 
state till the log-energy scale F ^ F^ = 1/S. Beyond 
this scale, the remaining odd bonds rapidly become much 
weaker relative to the remaining even bonds; the distri- 
bution of the even log-couplings /"=(C|r) = / dlP^iC, l\r) 
approaches some limiting distribution with a finite but 
large width, while distribution of the odd log-couplings 
-P°(C|r) = / dlP°{C,l\r) grows infinitely broad. In the 
RG language, the system rcnormalizes to some point on 
a line of RD fixed points (from this point of view, the 
RS states at S = represent critical points separating 
RD fixed points with opposite dimerization — see Fig. H). 
The corresponding joint distributions of the log-couplings 
and the lengths have been worked out in Ref. ^; here we 
only note that P'=(C|r) = ToiT)e-'^''°'^^'> with to(F) « 2S, 
while F°(C|F) = uo(F)e-«"«(r) with uo{T) « 2Se-^^^ . 
The ground state again consists of singlet pairs made up 
of one spin on an even site i and a second spin on some 
odd site j. Note however, that while i > j and i < j are 
equally probable in the RS state, in the RD phase with 
(5 > one almost always has j > i (with the exception of 
a few high-energy pairs of small spatial extent). 



B. Dynamic structure factor 

In this section, we summarize our calculations of the 
dynamic structure factor in the different regions of the 
phase diagram of spin-1/2 XXZ chains. Our approach 
has already been reviewed in general terms in Sec. [ID , 
and our calculations here represent one of the simplest 
examples of this approach at work. 

We begin by considering only the leading term in the 
perturbative expansion for the rcnormalizcd spin opera- 
tors; the results obtained in this manner give the correct 
leading behavior at low frequencies (some justification of 



this is given in Sec. Ill E| , where we discuss the role of 
higher-order terms). 



1. Random Singlet states 

The leading-order 'operator renormalizations' needed 
are particularly simple: the spin operator s remains un- 
changed for each of the 'surviving' spins and is effectively 
zero for each of the 'decimated' spins (i.e., spins that are 
already locked into singlets with other spins). 

Consider first S^^{k,uj); in our formulation of the sin- 



glet RG, Sec. Ill A , the following analysis applies to a 
general XXZ singlet state (i.e., remains valid so long as 
the ground state does not have lAF order). Consider 
two spins L and R connected by a strong bond {J-^,J^) 
in the renormalized theory with cutoff Ognai- The spin 
operators s£^^ connect the singlet ground state of this 
pair only to the triplet state \to) (with niz = 0), which 
is separated from the singlet state by a gap J-^. There- 
fore, the energy scale rifinai at which we stop the RG is 
^^finai — in this case (remember that the cutoff was 
defined as = max{J^}). We thus consider the renor- 
malized spectral sum 



'Sito-Er,,), (18) 



where the tildes remind us of the fact that this spec- 
tral sum now refers to the new Hamiltonian with en- 
ergy cutoff rifinai = w; this renormalized Hamiltonian 
has nr^ spins per unit length with the distribution of 
couplings and bond lengths characteristic of the fixed 
point to which the system flows in the low-energy limit. 
The sum Eq. ( [l8| ) is dominated by the excitations to the 
triplet state jto) of pairs of spins connected by the (renor- 
malized) bonds with J-^ ~ uj; these pairs are precisely the 
ones that are being eliminated at this energy scale. The 
corresponding matrix element for each such pair is simply 
(1 — e''''')/2, where I is the length of the bond connecting 
the pair; this allows us to write 



S'^^(/c,t^) 



7i(F^) J d/dC|l-e^"|2p(C,;|F^)(5(o 



(19) 



for w <C r^o in any RS state. 

The calculation of (fc, lu) is more involved since the 
gap to the relevant triplet excited state |ti) (with = 1) 
of a pair of spins connected by a strong bond (J-*-, J^) 
is now {J-^ + J^)/2. We consider each of the three cases 
(XX, XXX, and XXZC) separately: (l)ln the XXX case, 
the Heisenberg symmetry of the problem guarantees that 
gxx _ gyy _ gzz^ ^2) When the system approaches the 
XX point at low energies, we have <C J-*- implying 
that the relevant gap is approximately 12. Thus, to 
calculate (fc,ti^) we now have to stop the RG at the 
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scale rifinai = 2aj. ^From the Eq. (|T9|), it is clear that 
this leaves our answer unchanged except for the values 
of various non-universal scale factors. (3) The XXZC 
critical point needs special attention. In this case J^' j 
can have a range of values. As a result, the excited states 
that dominate the spectral sum Eq. (||) are not simply 
obtained by stopping the RG at any particular rjfinai and 
looking at the singlets forming only at this scale. Instead, 
for any fifinai G (0, Scj) there will be some singlets formed 
at this scale that will contribute to the spectral sum — 
namely, the pairs coupled by strong bonds with = 
f^finai and = 2lj — ilfinai- Notc that there is no double- 
counting here since we are considering only the pairs that 
are being eliminated at each energy scale. Thus, we have 

S+-{k,uj)r^ JdTdldDn{r)\l-e'''^\^P{0,D,l\T) x 

xS{uj-QQe-^{l + e°)/2). (20) 

Rewriting this in terms of the scaling probability dis- 
tribution V2 and using the delta function to do the F 
integral gives us 



S+-{k,uj) 



1 



wF„ 



XV2 



dl dD 



D 



(21) 



where we have defined D = D/T^, I ^ l/Tl, fc = fcF^, 
and 



ln(l + e^^^)-ln2 

f3 



(22) 



Now, since 1/' < 1; it is permissible to take the T^^ ^ 00 
limit of T^{D) before doing the D integral — in other 
words, we can replace T by 1 in the low-energy limit. 
The D integral can then be done trivially, and the final 
expression is identical in form to Eq. (E9h. More physi- 
cally, a given bond (J-*-, J^) is described fairly well (on 
a logarithmic scale) by one of these two couplings — we 
chose the characteristic scale to be J-^ . Now, the random 
anisotropy leads to an uncertainty | ln(J^/J-'-)| in 
the corresponding log-energy scale. This uncertainty is 
much smaller than the already existing typical spread in 
the log-energies or the typical log-energies themselves, 
which are both of order F. The leading behavior at low 
frequencies is therefore not affected. 

Thus, in the limit of low frequencies both S^^{k,uj) 
and {k,Lu) can be expressed in terms of the scaling 
probability distribution Vi as 

S{k,u;)r^^^ Lr|l„e'fcf 7,^(0,0. (23) 



Let us first focus on the regime \q\ = |fc — n/a\ ^ a^^. 
Note that all unsealed lengths I are odd multiples of the 



unit length a, and therefore e 



ikl 



The integral 



in Eq. ( pSf ) can be evaluated using the characterization 
of the function Vi{0,y) available in Ref. ^; the result is 
the following rather unusual scaling form at all three RS 
fixed points: 



o{k = \- q,uj) = 3 

a lyUjln {Qq/uj) 



(24) 



Note that we have suppressed the component labels on 
5(fc,w) as the two independent components obey the 
same scaling form, but with different values in general 
of the numerical constant A and the microscopic length 
scale ly. The universal function <I>(x) can be written as 



$(x) = 1 



cos(x) sinh(a;) + sin(a;) cosh(x) 
cos^(a;) sinh^(a;) + sin^(a;) cosh^(x) 



(25) 



The resulting S{k,uj) is shown on Fig. g. There is a 
fairly straightforward interpretation of the main features 
of this lineshape: The peak at g = (i.e., at fc = w/a) re- 
flects the predominantly antiferromagnetic character of 
the low-energy fluctuations; in our language, this is a 
direct consequence of the fact that the (renormalized) 
bonds all have odd lengths in units of a. The strongly 
damped oscillations with the period and the decay scale 
both of order F~^ express the properties of the distribu- 
tion of lengths of the strong bonds: both the average and 
the RMS fluctuation of this distribution of lengths are of 
order Ff,. 

While this result is interesting, one needs to analyze 
the effects of higher-order terms in the operator renor- 
malizations before accepting its consequences for possi- 
ble neutr on scattering experiments. We will argue in 
Sec. [II E that higher-order corrections do not modify the 
functional form Eq. (^5|) of the features in S{k, to) at fixed 
Lu but only add an "incoherent" background (of strength 
comparable to that of the features) and suppress the am- 
plitude of the features by a non-universal multiplicative 
factor of order one. 

A similar scaling function can be derived for the regime 
|fc| <C a~^. Repeating the above analysis gives 



S{k,Lu) = 



A' 



(|fc;,|i/2ln(17o/w)) , (26) 



with ^{x) — 2—^{x), and A' an order one numerical con- 
stant. This scaling function vanishes for fc ^ 0; for small 
k we have S{k,uj) ~ lyk'^ln{no/uj)/uj. We must there- 
fore consider the possibility that higher-order corrections 
may overwhelm this scaling result and render it irrele- 
vant. This is indeed expected to happen for {k,uj) 
away from the XXX point. However, we expect the scal- 
ing result to be valid quite generally for S^^{k,uj) — spin 
conservation guarantees that the higher-order corrections 
to S^^{k,uj) must also vanish as fe ^ (see Sec. IIIF for 
a detailed discussion of this point). 
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2. Dynamic structure factor in Random Dimer phases: 
'sharpness ' of the Griffiths regions 

Next, we consider spin dynamic structure factor in 
the XX and XXX Random Dimer phases introduced in 



Sec. [II A 3. The same approach as for the RS states goes 



over unchanged, and we write 

S{k,u})^n{r^) y"d/dC|l-e*"px 

X [P°(C, l\r^) + P^(C, l\r^)]S{cu - u:e-<) (27) 

for both S'^^(/c,w) and 5+~(fc,w) in both the XX and 
XXX RD phases (we are again being sloppy about the 
distinction between and as this can be absorbed 

in the definition of the non-universal scale factors that 
enter our expressions). 

Using the results of Ref. |^, it is a simple matter to 
obtain the full crossover from the RS-like behavior of 
the structure factor in the regime 1 ^ F^^ < Fa to the 
behavior characteristic of the RD phase in the regime 
F^ ^ F^. Here, we focus on the behavior in the regime 
F[^ S> F5, as this exhibits some rather unusual features. 
At these low energies, the even bonds dominate over the 
odd bonds, and the contribution of the odd bonds to the 
sum Eq. (g^) is negligible (we are assuming (5 > for 
concreteness). For wavevectors in the vicinity of k — n/a 
with |g| = |fc — 7r/a| <C S'^/ly (i.e., probing lengths larger 
than the correlation length ^av ~ Iv/S'^) we obtain 



c\6\^n- 



l/^RD 



l + cos(?„gF^/|5|)e-^''''r./l*'l' 

(28) 



where C and c are some order one constants, and we have 
chosen to write the power-law prefactor in terms of the 
dynamical exponent z^d (as far as our RG calculations 
are concerned = 2\S\ for small \S\ — however, the ef- 
fective value of 5 that enters this expression is expected 
to acquire a non-universal multiplicative renormalization 
from the high-energy physics, and the only reliable state- 
ment we can make is that z^^ ^ \^\ for small enough \S\). 
This result has a striking oscillatory structure (see Fig. |j) 
that is not suppressed significantly by the exponential 
factor, since \/F^/|(5|3 ^ T^/lf^l in the regime under 
consideration. This is best understood as a novel signa- 
ture of the sharply defined geometry of the rare Griffiths 
regions that contribute to the scattering at a given low 
energy (i.e., that are filtered out by their energy). More 
precisely, the average length of such regions is of order 
Iv^si^Lo/^s) = lv^uj/\S\, while the RMS fiuctuations in 
the length are only of order lyTg y/T^/Ts = Iv \/T^/\6\'^. 
Our results thus suggest that low-energy INS experiments 
would be able to pick up the sharply defined geometry of 
such Griffiths regions in the RD phases in one dimension. 

This feature of the Griffiths regions in one dimension 
was noted in Ref. ^ Sec. IVB in the context of the RT- 



FIM, where it was conjectured also that other proper- 
ties of such low-energy regions are likewise sharply de- 
fined: for example, in the disordered phase of the RT- 
FIM, the magnetic moment of the Griffiths regions with 
a given characteristic energy is sharply defincdcil and pro- 
portional to the (sharply defined) length of such regions. 
[In fact, similar "sharpness" is expected to hold for any 
bond "property" that "rides" on top of the singlet RG 
via recursion relation x = xi + + Tx2 when bond J2 is 
eliminated.] We expect to see a signature of this sharp- 
ness of the Griffiths regions also in the dynamic structure 
factor S^^ in the lAF Griffiths phase (see below) and 
also in the Griffiths phases of the one-dimensional RT- 
FIM (Sec. |V C| ). Finally, an interesting question, which 
we leave unanswered for now, is whether similar sharp- 
ness in the properties of the Griffiths regions at a given 
energy occurs and has observable consequences in higher 
dimensions as well, e.g., in the disordered phase of the 
d > 1 RTFIM. 



3. lAF Griffiths phase 



[II A 2 



in the lAF Griffiths 
the dominant low- 
are classical domain 



Let us first consider S^^{ k, m 

phase. As discussed in Sec ^ 

energy excitations in this phase 
walls. However, it is clear that such excitations do not 
contribute at all to S^^{k,uj), since they cannot be ex- 
cited from the ground state by the action of the operators 
like §1 , which conserve the total s^^^ . The leading excita- 
tions that do contribute to S^^ can clearly be identified 
in the RG picture with the rriz = excited states of pairs 
of super-spins, with each pair connected by a bond with 
J-L ijj and forming a singlet (note that this is true re- 
gardless of the value of the corresponding J^). Now, it 
is easy to generate a weak coupling of order cj in the 
lAF phase, since any typical region of length L will have 
an effective of order ^Qe~^^^ (and an effective typ- 
ically much stronger). What is more difficult is to isolate 
such a region from becoming a part of a larger cluster — 
otherwise this region can not support spin fluctuations at 
frequency lo. For this, we need two rare RS-likc segments 
(domain walls) with ^ w, one on each side of our (typ- 
ical) region. Thus, we need two domain walls, which are 
usually separated by a large distance of order to 
occur close to each other — the "density" of such occur- 
rences is '-^ (^2/ziAF^ 'pj-^g separation of the two domain 
walls — the length of the lAF-ordered cluster that they 
isolate — must be of order |lnijj|. More precisely, if the 
lAF-ordered cluster has length L, it can be thought of as 
consisting of the n ~ L 1^%^^^ strongly Ising coupled spins 
that are active at the crossover scale; the effective bonds 
connecting these spins at the crossover scale typically sat- 
isfy ln( J-*-/ J^) ~ — F^j^p. The requirement that the spin- 
flip coupling for this cluster is w flxes the length of this 
cluster to be i = ^tiF^^F^j^p, while the uncertainty in this 

length can only be of order /„i/F, 



^ OIAF 
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We are now ready to calculate S''^{k = n/a + q,LLj) 
in the regime \q\~^ ^ /^F^^^^, in addition to uj ^Siaf- 
The leading-order renormalization of the s| in the cluster 
RG is simple: s| is renormalized to (— l)-'s^ for each spin 
j that is active in some cluster c, and renormalizes to 
zero for every spin that forms a singlet. Assuming that 
such clusters "look" fairly uniform on the length scales 
larger than Iv^Sjj^-p: and adding up the contributions from 
all such isolated clusters with effective spin fluctuation 
frequency uj, we obtain 

-2/2;iAF 



1 



qHluj 
cos{qUT^/\5iA_-p 



1-2/21AF 



.r^/l^iAFl' 



(29) 



where C and c are some order one constants and the 
power of the <5iaf that appears in the prefactor has been 
fixed by demanding consistency with the off-critical scal- 
ing form 

5-(fc = - + = -4j* WA'"^^ 

a i^wli VJ^^iAF 

(30) 

with *(0,y) $(?/)• Note also that the overaU \/q^ 
dependence is a consequence of the fact that the spins 
contributing to the scattering have been taken to be dis- 
tributed uniformly over a sharply defined region (the 
cluster); we expect this to cross over to a much faster 
decay at large momenta (such that ~ I 
outside the range of validity of our scaling picture 

uite different for {k^uj). As we 



rLp) -ell 



IllC, the renormalization of the s 



± 



The situation is q 
shall see in Section 
spin operators is quite non-trivial, and we are unable to 
make an equally detailed prediction for 5"' . However, 
we expect that the matrix element for producing domain 
wall excitations with energies of order lo by the action of 
on the ground state is strongly suppressed as some 
power of cj, giving rise to a correspondingly small value 
for {k,oj) at small lo. 



C. Average local autocorrelations 

The same approach can be used to calculate average 
autocorrelation functions, and this section is devoted to 
a brief account of our results. 

We consider the local dynamical susceptibilities 



Xn (^) 



\{m\s^m^5{u-E^), 



(31) 



where a ^ z ov a — x. A knowledge of the low- 
frequency behavior of these susceptibilities can immedi- 
ately be translated into information about the long-time 
limit of the corresponding imaginary-time autocorrela- 
tion functions 



(32) 



1. RS states and RD phases 

As long as one is only interested in averages of such 
local quantities (over different realizations of disorder), 
it again suffices to consider only the leading-order spin 
operator renormalizations. We thus already have all the 
ingredients needed to calculate these average dynamical 
susceptibilities: our basic approach is familiar enough 
by now, and the relevant results of Ref. || for the renor- 
malized bond distributions have already been reviewed 
in Sec. Ill A. Below, we will be correspondingly brief. 



We first give our results for the average local dynamical 
susceptibilities and then translate these to results for the 
long-time behavioppf the corresponding average autocor- 
relation functions.E3 The leading behavior is the same for 
both a — z and a — x, so we drop all superscripts. 
For a bulk-spin, we obtain 



[Xloclav (^) 



^(Po^(Fj + P°(Fj). (33) 



For the critical RS states (P° — P°) we find 



[Xloclav (^) 
[Clod, (r) 



1 



lnw|3 



1 



llnrP ' 



(34) 



while off-critical — in the RD phases — we find 

\r 1 M 1^1' 



(35) 



Similarly, for an end-spin si of a semi-infinite chain 
(with j > 1) we obtain 



[Xllav ('^) 



Po^(r^)Po°(F^) 



(36) 



For the RS states we find 

[Xllav i^) ' 

and in the RD phases 



1 



u;| Inwp ' 
1 



llnrl 



(37) 



(38) 
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2. lAF Griffiths phase 

In the lAF phase, unhke in the singlet states, we need 
to make a distinction between and x^^. Consider 
first [xfcTclav ('^)- iFrom our previous discussion of the 
lAF phase, it is clear that, in the regime cu ^ ^^Siafj 
the dominant contributions come from lAF-ordered clus- 
ters of lengths ^ Ti^F^j^p (i-e., with effective spin- flip 
couplings of order w) that are isolated from the rest of 
the system by domain walls with ^ w on either side. 
^From the scaling picture of the phase, we get 

iXlocJav i^) ~ '^lAF ^l_2/^iAF ' 

l^^XAr)-Stl,j0^. (39) 

The analysis is more complicated for [xfocl i^)' ^^'^ 
we can only make a plausible estimate for this quan- 
tity. This is because the x and y components of the 
spin operators renormalize in a non-trivial way under the 
cluster RG. The origin of this difficulty may be seen as 
follows: Consider, for example, combining three spins 
S2, S3, and S4, connected by strong J| and J|, into a 
super-spin 5(234) ■ To zeroth order, all three operators sj, 

, and s^, renormalize to zero. To first order, and 
S4 renormalize to (sj)''^ — —s^J^/Ji ~ ^(234)-^"^3^/'^2 > 
(54)"'^ = — s^J^/J| — S(^34^2J^/J|, while Sg renormal- 
izes to zero to this order. Roughly speaking, the original 
spin-flip operators of the (active) spins have projections 
onto the remaining effective cluster spin-flip operators 
with components given by the ratio of the correspond- 
ing effective spin-flip couplings to the original spin-flip 
couplings. 

Now, the dominant contributions to [Xk)c]av i^) come 
from the low-energy (of order to) domain wall excitations, 
which are represented in the RG picture by the bonds 
with ^ Lo connecting the effective spins (clusters) in 
the effective theory with the renormalized cutoff lo. The 
matrix element for producing such an excitation by a 
bare spin-flip operator of a spin active in one of these 
clusters will be of order the corresponding , while the 
number of such spins contributing will be of order some 
effective "moment" Hx of this cluster. Because of the 
matrix element proportional to J"*", there will be a sig- 
nificant contribution only if this is also of order uo. 
As we have already seen, this can happen only if such 
an lAF-ordered cluster has length of order T^^ and is iso- 
lated from the rest by RS-like regions (domain walls) with 
^ on either side. We already know how to estimate 
the number density of such Griffiths regions. As far as 
the effective moment of such an lAF cluster is con- 
cerned, we can only make a crude estimate that bounds 
it from above by the number of spins that are active in 
this cluster: ^ixioj) ^ F^; however, we are unable to ob- 
tain the precise power of the logarithm that enters the 
energy dependence of the effective moment. We therefore 



leave out the logarithmic correction, and only write the 
dominant power-law part of our estimate: 

[Cloclav (-) - ^ • (40) 

D. Spin transport 

This section is devoted to a discussion of the dynami- 
cal spin conductivity cr' {uj) in the spin-1/2 XXZ chains. 
Our task here is to evaluate the Kubo formula Eq. (H) 
in the low-frequency limit. For the RS and RD states, 
we will use information available from the scaling solu- 
tions to the singlet RG recursion relations to achieve this, 
while in the lAF phase, we will use the scaling picture 
of the Griffiths phase we have developed earlier. Our re- 
sults for the dynamical conductivity are summarized in 
Figs. |l| and ^. 

1. Random singlet states 

We first need to work out the rules that govern the 
renormalizations of the current operators. Assume once 
again that is the strongest bond. We wish to work 
out perturbatively the renormalized operators tj^/2/3 that 
we trade in Ti/2/3 for, when we freeze spins 2 and 3 
in their singlet ground state (the other current opera- 
tors to the left and right of this segment are left un- 
changed to leading order by the renormalization). Now, 
note that these other operators have overall scale fac- 
tors in them that are nothing but the corresponding J-^ 
couplings. In order to be consistent, we clearly need to 
work out fi/2/3 correct to O(Jj^) (where J14 is the ef- 
fective bond connecting spins 1 and 4 after we freeze out 
spins 2 and 3) by adding the effects of virtual fluctu- 
ations to the projections of Ti/2/3 into the singlet sub- 
space. An explicit calculation gives the simple result 
that all three operators renormalize to the same oper- 
ator fi/3 = 72 = iJf4(s^S4 — s^s^)/2, which wc will 
denote henceforth by fi for consistency of notation. 

As we carry out the RG, the above result implies that 
the total current operator X]^=i entering Eq. (||) renor- 
malizes to '^jljfj, where j now labels the remaining sites 

of the renormalized system, and the Ij are the lengths of 
the corresponding renormalized bonds. [Note that this 
result makes sense physically and is a consequence of 
spin conservation: when a magnetic field with a uniform 
gradient is applied along the length of the chain, the ef- 
fective lengths Ij measure the "phase" along the chain of 
this "driving potential".] Consider two spins connected 
by a strong bond (J"*", J^) in the renormalized theory 
with cutoff rifinai- Siucc the current operator living on 
this bond connects the singlet ground state of the pair 
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only to the triplet state \to) separated from the singlet by 
a gap J"*", we choose ilfinai = w and consider the renor- 
malized spectral sum 

'^'M = ^El(m|D>.|0)|'<S(c.-^™). (41) 

rn j 

This spectral sum is dominated by precisely the \to) 
triplet excitations of pairs of spins that are connected 
by the (effective) bonds with J-^ — cu and are being 
eliminated at this energy scale; the corresponding ma- 
trix element is just loj/2, where I is the length of the 
bond connecting the pair. In the thermodynamic limit, 
we thus have 

a^uj)^"^^^ fdldCuj^f PiCJ\r^)S{LO - uje-'^) . (42) 
J 

This immediately yields our central result 

a'{io)^ICnslvH^o/Lo), (43) 

valid for a; <C f^o- Here, /Crs is an order one numerical 
constant, ly is the microscopic length scale defined earlier, 
and Og is the microscopic energy cutoff. Notice that this 
analysis holds equally well at all three RS fixed points, 
which differ only in the corresponding values of the non- 
universal scale factors. 

A brief digression is in order, before we go on to discuss 
this result: The real part of the dynamical conductivity 
can be related (on general grounds) to the behavior of 
the dynamic structure factor S^^{k,uj) near fc = 

<j'{io)^u:l^S^^k,uj); (44) 

this can be checked by comparing directly the corre- 
sponding spectral sums and noticing that the action of 
the two operators T = J^j and V = J2j J'^j on the 
eigenstates of the Hamiltonian Ti are related through 
T = i\H, V]. It is easy to check, using the scaling form 
Eq. (p6|), that our result for the conductivity is consis- 
tent, as it must be, with our previously derived result for 
the dynamic structure factor. 

Going back to Eq. (^), we see that cr'(w) diverges log- 
arithmically for small w in the unusual 'spin-metal' phase 
controlled by the XX fixed point as well as at the critical 
points (XXX and XXZC) separating this phase from the 
'insulating' phase with Ising antiferromagnetic order in 
the ground state. Note that this 'metal-insulator' tran- 
sition has the curious feature that the quantum critical 
points separating the conducting phase from the insulat- 
ing phase have the same T = transport properties as 
the conducting phase. 

2. lAF Griffiths phase 

On the insulating side, we expect a'{uj) to be sup- 
pressed below the crossover scale ^2^;^^; the dominant 



contributions for lo ^^iSiaf come from some rare regions 
that contain long finite segments locally in the 'metallic' 
phase. 

We begin by providing a rough estimate of these con- 
tributions to cr'(ti'): In our sample, consider a (large) 
region of length L locally in the RS phase; the number 
density of such regions is roughly , with some p < 1 
(which depends on the distance from the transition). If 
these regions are effectively isolated from the rest of the 
system, the average power absorption per spin in each 
such region is proportional to the finite-size conductivity 
calculated in Appendix |^ 

W = icrRs(w, L) ~ exp(-c| lncj| V^) , (45) 

where we have assumed that L, although large, satisfies 
L ^ jlnwp (this assumption will turn out to be self- 
consistent). The total power absorbed in the sample is 
then obtained by summing over all such regions: 

a'{uj)r^ JdLp^ L^/^exp{~c\\nuj\^/L). (46) 

Evaluating this integral by a saddle-point method, we 
find that the lengths that dominate are of order | In a; | 
(our assumption about the lengths is thus valid), and 
arrive at the following estimate 

a'{Lj) r^LU°'\\nLj\^ , (47) 

where a = a((5iAF) > is a continuously varying expo- 
nent vanishing at the transition. While this argument 
is suggestive, we find it more convincingcS to take an 
alternative route based on the scaling picture we have 
developed earlier for the lAF phase — this has the added 
advantage that it allows us to relate the exponent a to 
the dynamical exponent ^(Jiaf)- This is what we turn 
to next. 

We have already seen that the most numerous low- 
energy excitations in the lAF Griffiths phase are domain 
walls, with the integrated density of states ^ ^i/^iaf^ 
Such classical Ising excitations, however, do not con- 
tribute to the dynamical conductivity. The dominant 
contributions come from lAF-ordered clusters of lengths 
L ~ r^^r^j^p (i.e., with effective spin-flip couplings of 
order cu) that are isolated from the rest of the system 
by domain walls with ^ to. Remembering that the 
number density of such Griffiths regions is ^ ^2/ziaf^ 
and noting that the corresponding "phase lengths" are 
of order L ^ | lnti^|, we immediately obtain Eq. ( ^ ) with 
a — 2/ziAF- More formally, we sum over the possible 
separations of two such domain walls, with the constraint 
that the typical lAF-ordered region isolated by the two 
has significant spin fluctuations at the characteristic fre- 
quency Lo: 

a'(cj) - ^ / dLuj^L^6{iJ ~ noe-^-^) . (48) 
^ J 

We thus obtain for the dynamical conductivity in the lAF 
phase 
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= JCjAFlvi^/ftoP^'^^ ln\rto/Lo) , (49) 

where /Ciaf is a numerical prefactor that depends contin- 
uously on (5iAF- The scaling of /Ciaf with 6iaf for small 
(5iAF can be obtained by demanding consistency with the 
ofF-critical scaling form for the conductivity 

cr'(w) = K-Rslv ln(f7o/w)i;iAF(ri^/r5iAF) , (50) 

which immediately implies that /Ciaf ~ '^iaf'''^^^ ^ ■^iaf- 

3. Dynamical conductivity in RD phases 

We now calculate the dynamical spin conductivity in 
the XX and XXX Random Dimer phases. Here, the 
same singlet RG can be employed all the way across the 
crossover scale Tg = l/\S\, and into the energy regime of 
a well-developed RD phase. The dynamical conductivity 
is given by the same expression Eq. (^) as for the RS 
states: we simply add contributions from the even {P") 
and the odd (P°) bonds in complete analogy with the 
calculation of the dynamic structure factor. Using the 
scaling solutions of Ref. ^, it is quite simple to calculate 
the full scaling function for the dynamical conductivity 

a' {to, S) = /Crs/„ ln(r!o/^) Srd {\6\ ln(r!o/w)) ■ (51) 

Here, we restrict ourselves to noting that I]RD(a^) ^ 
const for X <C 1, while for a; 3> 1, SRD(a;) scales as 
^Rnix) xe~^^ . Thus, at frequencies uj well below the 
crossover scale fls, we have 

a'iio) = ICnnlv{Lo/no)'^'"° ln\no/to) , (52) 

with the numerical prefactor /Crd ~ |^| and the dynam- 
ical exponent zrd ^ |^|~^ for small \5\. We can now 
interpret this form directly in terms of the rare regions 
that dominate the conductivity: Assume, for concrete- 
ness, that (5 > 0, i.e., that the even bonds are dominat- 
ing; the main contribution to the dynamical conductivity 
at frequency uj ^ ftg then comes from the even bonds 
with effective — lo. Such weak even bonds are gen- 
erated only across rare long regions that are locally in 
the opposite dimerized phase, and these are precisely the 
regions that dominate the low-energy density of states 
and thus determine the dynamical exponent 2:rd(5) — 
this explains the factor uj^/^^^o jj-^ ( |52|) . Moreover, 
all such bonds have a well-defined length proportional to 
ln(r2o/i^), which explains the Iu^IQq/uj) in Eq. (|5^). 

4- Perspective: spinless interacting fermions with 
particle-hole symmetric disorder 

To put these transport properties in perspective, we 
recall that the spin-1/2 XXZ chain is equivalent, via 
the usual Jordan- Wigner transformation, to a system of 



spinless interacting fermions with particle-hole symmet- 
ric disorder. More specifically, we write the spin opera- 
tors = ± is^ in terms of fermion creation (annihi- 
lation) operators (cj) as 

4- 11(1- 2%-') c], 

j' <j 

n(l-2n,0c,, (53) 

j'<j 

while Sj = rij — 1/2 (here nj = CjCj is the fermion number 
operator at site j). In this language, Tixxz can be written 
as 

^ = E + + ^^■('^J- - V2)(ri,+i - 1/2) 

(54) 

with tj = Jj~/2 and Vj — Jj. The coupling thus con- 
trols the strength of the nearest-neighbor particle-hole- 
symmetric repulsive interaction between the fermions. 
The lAF phase that obtains for large corresponds to 
a charge density wave state stabilized by interactions. 
In the absence of interactions (XX chain) we obtain a 
free-fermion random-hopping problem at zero chemical 
potential. 

This free-fermion problem has been extensively stud- 
ied in the past, and is known to have rather unusual lo- 
calization properti| eS || d ue to the additional particle-hole 
symmetry present.EJcJ For instance, an elementary cal- 
culation immediately reveals that the zero-temperature 
average Landauer conductance [f^ijav of a finite segment 
of length L connected to perfect leads scales as [^ijav ~ 
in sharp contrast to the usual exponentially- 
localized behavior in one dimension; the corresponding 
conductivity, of course, scales as y/L. Now, the strong- 
disorder RG predicts that lengths scale as the square of 
the logarithm of the energy scale in the low-energy effec- 
tive theory describing the XX-RS state — our result for 
the dynamical conductivity is thus consistent with the 
elementary Landauer calculation (see also out, explicit 
finite-size scaling calculations of Appendix |^c3 Notice, 
however, that our approach is not limited to the non- 
interacting case. It allows us to reliably treat the effects 
of interactions, and follow the dynamical conductivity 
through a 'metal-insulator' transition that is driven by 
strong interactions in the presence of strong disorder. 

5. Numerical study of the dynamical conductivity at the XX 
fixed point 

At the XX point, the Hamiltonian Eq. ( [5^ ) describes 
non-interacting fermions with random hopping ampli- 
tudes, and we are essentially faced with the problem 
of finding the low-energy eigenvalues and eigenstates of 
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the corresponding single-particle Hamiltonian (an L x L 

matrix operator) H = Ylpl tjilJ + 1)01 + b'Xi + 
which defines the Schroedinger equation for this prob- 
lem. Any fermionic state can then be represented as 
a Slater determinant of the corresponding (normalized) 
single-particle eigenstates |<^^) with eigen-energies e^. In 
the single-particle language, the Kubo formula for the 
conductivity (t'{lu) at zero chemical potential and at a 
finite temperature T reads 



ujL 

X [/(e^i) - /(eM2)]'5(w - £^2 + e^i) ' (55) 

where T{j) = itj{\j){j + 1\ ~ \j + is the current 

operator on the link j and /(e) = l/(e'^/"^ -I- 1). [This 
version of the Kubo formula will also prove useful when 
we analyze the full temperature dependence of the dy- 
namical conductivity in Sec. VI.] 

Here, we test the T ~ predictions by evaluating 
ct'(w) using exact numerical diagonalization of finite sys- 
tems. The results of such calculations for system sizes 
L = 128, 256 and 512 with the hopping amplitudes tj 
drawn independently from a uniform distribution over 
[0, 1] are shown on Fig. ^, where we have averaged over 
100, ODD samples for each L. In an infinite sample we ex- 
pect the conductivity to diverge logarithmically, but with 
the system sizes studied here, we cannot quite probe this 
infinite-sample regime 1 <C \n(rio/uj) ^ \/L — rather, we 
are in the regime 1 <C ln(rio/t^) ^ VL. Nevertheless, the 
numerical results of Fig. ^ clearly show that the dynami- 
cal conductivity is increasing as the frequency is lowered 
all the way to the crossover scale \ii{flo/uj) ~ a/L, thus 
supporting our claim that (t'(w) diverges logarithmically 
at low frequencies. 

For a more detailed test of our theoretical results, we 
need to quantitatively analyze the effects of a finite sys- 
tem size on our predictions for the dynamical conduc- 
tivity. The calculation is summarized in Appendix 
Here, we only note that this analysis allows us to write 
the following scaling form for the conductivity 



cr'(w,i) = ?„ln(f7o/w)e (/„ \n{no/Lj)/L) 



(56) 



the scaling function Q is characterized in Appendix 
and the above result is expected to hold for large 
enough L and ln(f2o/<^) (with no restrictions on the ra- 
tio In^ / 1-^) / L) ■ However, the numerical results cannot 
be compared directly with this scaling result since it as- 
sumes that the distribution of bond lengths has reached 
the form characteristic of the XX fixed point, which is 
not the case for the sizes that we can diagonalize nu- 
merically: the "length part" of the distribution P(^, ZjF) 
is still evolving towards the corresponding scaling form 
from the initial condition P(C, = e^'^S{l — 1). Nev- 
ertheless, we can compare the results of the exact numer- 
ical diagonalization with the (formal) predictions of the 
RG for the same systems. This can be done by either 
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FIG. 6. Scaling plot of the dynamical conductivity cr'(cj, L) 
at the XX point for finite systems of sizes L — 128, 
256, and 512, with free boundary conditions, calculated by 
exact numerical diagonalization and — for exactly the same 
systems — by the finite-size RG analysis of Appendix ^ 
Free-fermion hopping amplitudes tj are drawn independently 
from a uniform distribution over [0, 1], and the (bare) phase 
lengths are set to Ij =1. We used Qo ~ 2 corresponding to 
the initial energy cutoff in the equivalent spin system; also, 
we used Iv ~1 corresponding to the microscopic length scale 
in the problem. The agreement of the RG predictions with 
the results of the exact diagonalization is fairly good (given 
the not so strong initial disorder), and the dynamical con- 
ductivity is indeed increasing all the way to the (finite-size) 
cross-over scale ln(no/i^) ~ \/ L/lv Note that with the sizes 
studied we can only partially access the bulk scaling regime 
ln{Q,o/Lo) ^ y^L/lv, in which we expect a'{u,L) ^ ln(S7o/<jj) 
and which on the plot is towards the left of the horiz ontal axis. 
Also note that in the opposite regime \n{Q.(_)/u)) > L/l^, un- 
like in the bulk regime, a'{u>,L) is not self-averaging; in this 
regime, the plotted average over different samples represents 
roughly the distribution of the lowest gap in the system. Inset 
shows how the L^oo scaling form (thick line) is approached 
by the finite-size o-'{uj, L) calculated from the RG for the given 
initial conditions (the lines plotted here are the same as in the 
main panel). Note that the vertical scale is set by the (nu- 
merical) prefactor of the scaling function in the bulk scaling 
regime (where cr'(Lj) ~ -^\n{Q,Q/ijj)) 



running the RG on the same samples or by evaluating 
the analytical (within the RG) expression, given by the 
inverse Laplace transform Eq. ( |A9| ) for these initial con- 
ditions. In Fig. ^, we compare the RG result obtained in 
this manner with the numerically evaluated conductiv- 
ity. Given that the initial disorder is not very strong, the 
agreement of the RG predictions with the (j'{u),L) from 
the exact diagonalization is fairly good. 
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E. On validity of results 

So far, our calculations have relied on the leading-order 
renormalizations of the spin operators; in this subsec- 
tion we will try to justify validity of this approxima- 
tion. We will not address the corresponding question for 
the RG itself — this has been analyzed with great care in 
Refs. |6|and^, and we have nothing to add here. Instead, 
we focus on issues specific to our calculation of dynam- 
ical quantities. Here, we provide a (partial) justifica- 
tion of our leading-order results by analyzing the effects 
of the first corrections to the leading-order expressions 
for the renormalized operators — this can be done consis- 
tently within the framework of the RG approach. Any 
consistent analysis of further corrections would require 
that we also consider higher-order corrections to the Hp:, 
rules themselves, and we stop well short of doing that.EII 

As an illustrative example, we consider the dynamic 
structure factor in the RS states. Our leading-order cal- 
culations used only the zcroth-order result for the renor- 
malized spin operators. The renormalized operators can 
also be easily worked out to first order — these were con- 
sidered in Ref. ^in the discussion of typical correlations. 
When a pair of spins 2 and 3 connected by a strong bond 
is frozen into a singlet state, the neighboring spin oper- 
ators 5*1 and 5*4 do not change even to first order, while 
the spin operators S2 and S3 renormalize to 

/ zNcff _ ( z\cS - "^l „z I J3 z 
[^2) — ~{^3) ~ ~2J^ 2J^ 

(4)^ff = -(4)-ff = jL,+ + _^4. (57) 

Ji + J 2 J2 +J2 

Thus, the decimated spins S2 and S3 obtain small "com- 
ponents" of order J1/J2, J3I J21 onto the neighboring 
spins si and S4. As we run the RG and renormalize 
down to scale F, the system consists of nr active spins per 
unit length, separated from each other by "dead" regions 
(with lengths of order F^) of decimated spins. Each dec- 
imated spin r in the dead region between two remaining 
active spins j and fc (where j and fc are nearest neigh- 
bors in the effective Hamiltonian at scale F) will have 
some components Crj and Crk on spins Sj and Sfc . ^From 
the point of view of calculating the spectral sum Eq. (|^), 
each active spin j acquires some (fc-dependent) effective 
moment flj{k) coming from all decimated spins with non- 
zero components on sj: 

fij{k) = l+ aje*^*''"^\ (58) 

ii <r<i2 

where the sum is over all previously decimated spins r 
between the effective neighbors ii and 12 of the spin j, 
ii < j < 12- The components Crj are simply the ground- 
state correlations (Si-Sj); such typical correlators decay 
as a stretched exponential — [InC^jJav ^ \r — Note 
that the characteristic length scale for this decay is the 



microscopic lerigth scale Z„. It is thus clear that the sum 
over r in Eq. (jSg) converges quickly, and the renormal- 
ization of the moment flj away from its bare value of 1 
comes mainly from the nearby spins that were decimated 
early in the RG. This renormalization is of order one, but 
only weakly fc-dependent. 

We now analyze the consequences of this renormal- 
ization of the moments for the two scaling forms of the 
dynamic structure factor derived earlier in the limit of 
low frequencies, one in the vicinity of fc = 7r/a, and 
the other in the vicinity of fc = 0. First, consider 
k — n/a + q, with \q\ <C ly^- For such small values 
of q, we can neglect the g-dependence of the moments 
and evaluate them a,t k — n / a. To evaluate the spectral 
sum Eq. (^), we need to add up the contributions com- 
ing from the strong bonds at scale rifinai- Each strong 
bond contributes |/iL -|- /2/je*^'p, where and are 
the moments (evaluated at fc = "'/'j) of the two spins 
connected by this strong bond. We can now proceed in 
two steps: First, we fix I and average over the moments 
of all strong bonds with a given length. This gives us a 
quantity ci -l- C2II -l- e*"^'!^ which we now need to average 
over the length distribution of the strong bonds; here ci 
and C2 are now some fixed numbers of order one, since 
we expect that the main renormalization of each moment 
comes from few nearby spins and is roughly independent 
of the lengths of the adjoining bonds. Thus, we see that 
Eq. (^), with $ given by Eq. (|25|), indeed describes the 
dynamic structure factor for fc close to n/a and fixed low 
uj; the higher-order corrections renormalize the overall 
amplitude by a factor of order one, and also produce an 
"incoherent" background of a comparable strength that 
depends only weakly on fc (i.e., that changes significantly 
only when fc is changed by an amount of order l^^). 

For k <C l^^ (i.e in the scaling regime near fc — 0), the 
discussion is very similar; each strong bond contributes 
I Mi — A_Re**''P, where the moments are now evaluated at 
fc = 0. This again gives us a quantity C1-I-C2 1 1-1-6**' p to be 
averaged over the length distribution of the strong bonds. 
In general (away from the XXX point), we now have 
to consider the S^^(k,uj) component separately from the 
(fc,w) component, since the total Sj^j conservation 
constrains the constant ci to be identically zero for the 
case of S^^{k,uj). Thus, in the case of S^^{k,uj), higher- 
order corrections only produce an order one renormaliza- 
tion to the overall scale of our scaling result Eq. (p6|); of 
course, there will be additional corrections, but these will 
vanish faster than the scaling result in the low-frequency 
limit. In the case of S'^'_(k,Lu), an inspection of the 
renormalization rules Eq. (|57| ) shows that to this order ci 
will be zero for (fc, uj) as well, even in the absence of 
full Heisenberg symmetry; however, this is not expected 
to be true in general (to all orders), and we expect a 
small but non-zero background to be present in the gen- 
eral case. Thus, (fc, u>) near fc = will in general 
consist of two parts: the scaling part given by Eq. ( |57| ) 
with an order one non-universal overall scale (this part 
vanishes as ~ fc^ for small fc), and a non-scaling weakly 
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fc-dependent additive background of the same order as 
the scaHng part. 

The above arguments typify the general logic behind 
our justification of the leading-order results for all of 
our calculations; in some cases such a programme can 
be carried out analytically (e.g., for the avera ge bound- 
ary spin autocorrelations — see also Sec. VD), while in 
other cases we have to be satisfied with arguments like 
the ones presented above. Such arguments can also be 
bolstered by numerically implementing the higher-order 
operator renormalizations to calculate corrections within 
the RG to our leading order results (indeed, we have 
confirmed that such a numerical check for S{k,u!) in the 
Heisenberg model is in qualitative agreement with the 
arguments presented above). 



IV. TRANSPORT IN STRONGLY RANDOM 
SPIN-1 CHAINS 

A. Singlet RG description of the phases: A review 

The strong-randomness quantum critical point, which 
controls the transition from the Haldane state to the 
Random Singlet state in the spin-1 chains, and the im- 
mediate vicinity of this critical point, can be analyzed 
by a somewhat extended RG procedure introduced in 
Refs. |lj and ||, or by a variant of the same used in Ref. ||. 

The basic idea is to replace the original spin-1 chain 
by an effective model that is argued to describe the low- 
energy physics of the original system — as we shall see 
later, this effective model can be made plausible by think- 
ing in terms of a bond-diluted chain (it is also possible 
to arrive at essentially the same model by starting with 
a random antiferromagnetic spin-1 chain and using the 
approximate^ RG procedure of Ref. H). This effective 
model is written entirely in terms of spin-1/2 degrees 
of freedom coupled by nearest-neighbor Heisenberg ex- 
change couplings. All even bonds are always antiferro- 
magnetic and are drawn from an appropriate distribution 
of positive bonds, while odd bonds can be of either sign 
and are drawn from a different distribution. 

This effective model can be analyzed using the exten- 
sion of the singlet RG introduced in Refs. |18| and |^. One 
begins by looking for the largest antiferromagnetic bond 
in the system, say J2 connecting spins 2 and 3; this de- 
fines our bare energy cutoff Qq. Further analysis can be 
split into three cases: (i) H the bonds adjacent to the 
largest AF bond are smaller in magnitude, the two spins 
are frozen into a singlet state and an effective coupling 
Ji4 is generated between spins 1 and 4 exactly as in the 
singlet RG for the spin-1/2 chain, (ii) If both adjacent 
bonds are larger in magnitude than J2, then spins 1 and 
2 and spins 3 and 4 are first combined to make effective 
spin-1 objects (since in this case Ji and J3 are necessarily 
ferromagnetic), and these effective spin-1 degrees of free- 
dom are then frozen into a singlet state, generating an 



effective coupling J05 = 4Jo J4/3J2 between spins and 5. 
(Hi) If only one of the adjacent bonds, say J3, is larger in 
magnitude than J2, then spins 3 and 4 are first combined 
into an effective spin-1 object. The system is then frozen 
into the subspace in which spin 2 and this effective spin- 
1 object are coupled together to form an effective spin- 
1/2 object which we label S2 for consistency of notation. 
The corresponding renormalized couplings are given as 
J12 ~ — Ji/3 and J25 = 2J4/3. This procedure is now it- 
erated with the energy cutoff fl being gradually reduced. 
It is important to note that there is no inconsistency in 
leaving untouched ferromagnetic bonds J < —Q that are 
not adjacent to any antiferromagnetic bonds at the cutoff 
scale — we could equally well have combined all pairs of 
spins connected by such strong ferromagnetic bonds into 
effective spin-1 objects at the cost of cluttering up our 
notation. 

A detailed analysis oLtliis iterative procedure can be 
summarized as follows :E3iJ'EI Let F = ln(rio/f^) and let 
nr be the fraction of active spins at log-cutoff F. For the 
even bonds, we introduce the distribution P(C|F) of the 
corresponding logarithmic couplings C = In(f2/J). For 
the odd bonds, let A^(F) be the fraction of odd bonds 
at scale F that are strongly ferromagnetic with J < — il; 
for large F, the remainder of the odd bonds are sym- 
metrically distributed about zero and are therefore de- 
scribed by a distribution for \J\ that we characterize 
by the distribution Q(C|F) of the corresponding loga- 
rithmic couplings C = l'^(^/l'^l)- When W, the width 
of the distribution of the log-exchanges in the original 
spin-1 Hamiltonian Eq. (|6|), exceeds a critical value Wc, 
the system is in a spin-1 Random Singlet phase. In the 
language of the spin-1/2 effective model, this RS phase 
is described by a fixed point with P(C|F) = F~^e^^/^, 
N{T) = 1, nr ~ 1/r^ and Q(C|r) = Qoe-Q°'^ for large 
F (Qo is some non-universal 0(1) number). As W is 
decreased, the system undergoes a quantum phase tran- 
sition to the so-called Gapless Haldane (GH) phase; both 
the quantum critical point and the GH phase in the vicin- 
ity of it are still controlled by strong-disorder fixed points. 
At the critical fixed point (which is an infinite-disorder 
fixed point) we have P(C|r) = Q(C|r) = 2T-^e-^'^/^, 
nr ^ l/F^, and A^(F) = 1/2. The GH phase in the 
vicinity of the quantum critical point is controlled by 
a line of fixed points; each point on this line is char- 
acterized by some constant Pq (which depends on the 
strength of disorder W). At a point labeled by Pq, we 
have P(C|F) = Pqc-^"'^, Q(C|F) = Qo(r)e-Q»(r)^ where 
Qo(r) ~ e--P°r, N{T) 0, and nr ~ Pge-^"^. The 
continuously varying Po(W) vanishes at the transition 
as Pq ~ {Wc — WY''^, where is the correlation length 
exponent obtained in Refs. |l8| and |^ — the GH phase 
is thus similar to the dimerized phases of the spin-1/2 
chains. 
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B. Spin transport 



2,. Dynamical conductivity 



1. Doing calculations in the effective model 

Before we calculate anything, we need to describe how 
we think about the spin transport in this case. This is 
somewhat non-trivial, for we are working in an effective 
model of spin-1/2 degrees of freedom, and some thought 
is required to decide what is the correct quantity to cal- 
culate. 

For this, we go back for a moment to the original ran- 
dom spin-1 chain and review an intuitive construction 
that leads to the effective model in terms of spin-1/2 
variables only. Consider the case of dilute randomness ,ta 
that is, consider a uniform spin-1 chain with a small frac- 
tion of very weak bonds that effectively break the chain 
into pure finite segments weakly coupled with each other. 
The low-energy effective degrees of freedom of such a seg- 
ment are two half-spins localized near the two edges of the 
segment — these are the spin-l/2's of the effective model. 
The coupling of the edge spins on neighboring segments 
is given roughly by the original coupling of the two seg- 
ments, and is always antiferromagnetic — these are the 
even bonds of the effective model. On the other hand, the 
coupling of the two edge spins of the same segment can 
be either antiferromagnetic or ferromagnetic depending 
on whether the length of the segment is even or odd — 
these coupling are represented by the odd bonds in the 
effective model. 

We now need to express dynamical properties of the 
system in terms of these effective spin-1/2 degrees of free- 
dom. In particular, we want to analyze the low- frequency 
power absorption when an oscillating magnetic field with 
a uniform gradient is applied to the system; this will give 
us the dynamical conductivity cr'(w). Since the mag- 
netic field couples to the conserved 'charge' in the sys- 
tem, the corresponding current operators that we need 
to use when working out the Kubo formula for the effec- 
tive model are uniquely determined by spin conservation: 
The current operator on the odd bonds connecting the 
edge half-spins si and S2 of the same segment (which rep- 
resents the total spin current operator of this segment) 
is T = J12I12S1 X 32] here J12 is the corresponding ef- 
fective coupling and I12 is some effective phase length 
that we expect to be given roughly by the length of the 
segment. Naturally, the current operators on the even 
bonds connecting the edge half-spins of the neighboring 
segments have a similar form (the argument in this case 
is even simpler: one only needs to know that the true 
edge spin-1 operator of a segment "projects" onto the 
corresponding effective edge spin-1/2 operator with an 
amplitude of order one). Note that the precise values of 
the phase lengths in the initial effective model (for the 
dilute spin-1 chain) are not important, since at still lower 
energies we expect the distributions of couplings and the 
corresponding bond lengths to approach some universal 
distributions characteristic of the appropriate fixed point. 



Having identified the appropriate current operators in 
the effective problem, we now work out the rules that 
govern their renormalization in the RG scheme used to 
analyze this effective model. As in the spin-1/2 case, and 
as discussed above, we write the part of the total current 
operator (in the spectral sum Eq. (||)) that is associated 
with a given bond + 1) in the form Ijfj where Tj is 
the usual bond operator fj — JjSj x sj+i and Ij is the 
appropriate phase length. We can then follow renormal- 
izations of the needed operators by keeping track of the 
phase lengths, in addition to the various bond-strengths. 
Unlike the spin-1/2 chains, these phase lengths need not 
equal the physical distances between the corresponding 
spins; in fact, even the physical position of an effective 
half-spin often cannot be specified unambiguously, as, for 
example, when this half-spin appears as an effective dou- 
blet formed by combining (via a strong AF bond) an 
effective spin-1 (which is an intermediate construction in 
the Hyman-Yang RG rules) and a neighboring spin-1/2. 
In such cases, our rules can actually be used to assign 
some meaning to the physical position of such an effec- 
tive half-spin. 

The rules for the phase lengths can be easily 
stated: In the cases (i) and (ii), when in the final 
step we form a singlet from either two spin-1/2 ob- 
jects or two spin-1 objects, the phase length of the 
new effective bond is simply the sum of the phase 
lengths of all the bonds that are eliminated. In the 
case (in) the phase lengths associated with effective 
bonds J12 and J25 are I12 = + (4/3)^2 + (2/3)/3 and 
I25 = k - (1/3)^2 + (1/3)^3 respectively. 

The rules for the phase lengths in the case (Hi) are 
somewhat unusual — for example, negative phase lengths 
can be produced. Note, however, that there are many 
factors that prevent this from happening too often, and 
the phase lengths will in many instances coincide with 
the corresponding geometrical lengths: decimations in 
the cases (i) and (ii) tend to "correct" deviations of the 
phase lengths from the geometrical lengths, and in both 
the RS and GH phases there are simply no decimations of 
type (Hi) at low enough energies. Also, the lengths I2 and 
I3 in the above rule for the case (Hi) are the lengths of the 
strong bonds that are eliminated and are therefore usu- 
ally smaller than the lengths li and I4 of the more typical 
bonds. Finally, one can argue generally that the phase 
positions of the spins as dictated by the phase lengths 
have to agree — at least roughly — with their geometrical 
positions as inferred from the order of the (remaining) 
spins in the chain (i.e., from the spin labels) O All of 
this implies that the phase lengths are roughly given by 
the geometrical distances between the spins; in particular 
their scaling with F is given by the inverse of the density 
of the remaining spins, / n{T)~^. 

We can now immediately deduce behavior of the dy- 
namical conductivity in the different phases exactly as 
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in our previous calculations for the spin-1/2 model; as 
the method remains the same, and the relevant details 
about the statistics of the fixed point Hamiltonians have 
already been summarized, we merely state our results. 

In the RS phase the same result Eq. ( |4^ ) applies, as 
is true for an RS state of an arbitrary-S* spin chain at 
strong enough randomness. 

At the critical point separating the RS phase from the 
GH phase, we find 

a'H =/CHY/.ln'(r!oM, (59) 

which is a stronger divergence than in the RS phase (note 
that this difference from the result in the RS states can 
be traced to the fact that the density of remaining spins 
behaves as nr ^ l/F^ at the critical point, in contrast to 
the 1/r^ decay of the corresponding quantity in the RS 
states) . 

Finally, In the GH phases parametrized by Pq we find 

a'iuj) = /CGH/.(c^/r!o)'/"'^" ln\no/Lu) , (60) 

where we have introduced the continuously varying dy- 
namical exponent zqh = Pq^, and /Cqh is an order one 
numerical prefactor which goes to a constant as ^ Wc 
(note that the factor ln^{^lo/u>) appears for exactly the 
same reasons as in the RD phases of the spin-1/2 chains: 
the lengths of the singlets that are decimated at scale u! 
are roughly ~ ln(rio/w))- 



A detailed analysis of this procedure was given in 
Ref. ^, of which a summary follows: Define the log- 
couplings Pi = \n{il/hi), C,i = \n{n/Ji), and also the 
log-cutoff r = \n{ilo/n); also, let np be the number den- 
sity per unit length of the (remaining) clusters at scale 
r. The essential feature of the RG near the critical point 
is that the distributions of the log-couplings R{(3\r) and 
P(C|r) become broader and broader as the energy cutoff 
is lowered; the RG flows are characterized by a special 
family of scahng solutions with R{f3\T) = i?o('r)e"-^°(^)'^ 
and P(C|r) = PoiT)e-^°^'^^'^. At the critical point, 6 = 0, 
we have Ro{T) = Po(T) = l/T; thus the widths of the 
two distributions grow without limit, and the number 
density decreases as np ^ l/F^. Also, magnetic moments 
of the clusters scale as /i ~ r-^, with <j)= {1 + V5)/2. In 
the disordered phase, (5 > 0, beyond the crossover scale 
= 1/\S\, the width of the field distribution saturates, 
with Ro {T) w 26 for T > T^, while the width of the bond 
distribution grows without limit, with Po(r) ~ 2Se^^^^ . 
In the ordered phase, 6 < 0, the situation is reversed: 
i?o(r) ~ 2|(5|e-2|'5|r and Po(r) « 2\S\ for T » Ts. In 
both phases, we have np ~ j^pe"^'''!'". Note also that 
the clusters that are being eliminated at scale F ^ 
all have a fairly well-defined length of order |(5|^^F and 
magnetic moment of order |i5|^~'^F. 

B. Average autocorrelations 



V. DYNAMICS IN THE RANDOM 
TRANSVERSE FIELD ISING CHAIN 

A. Strong-disorder RG description of the phases: A 



The strong-randomness cluster RG of Ref. ^, from 
which the low-energy long-distance behavior of a system 
near the critical point {\S\ <C 1) can be obtained, pro- 
ceeds as follows: One finds the largest coupling in the 
system, with energy fig = max{/ij, Jj}. If the largest 
coupling is a field, say /12 on spin 2, this spin is frozen 
into the af = -f 1 ground state of the local field term 
and is eliminated from the system leaving an effective 
coupling J13 = Ji2t/23/^2 between the neighboring spins 
1 and 3. If the largest coupling is an interaction, say 
J12 between spins 1 and 2, the two spins are combined 
into one new spin — a cluster — with an effective spin vari- 
able (7(12) [representing the two classical minimum energy 
states (jf = (7| = ±1] and an effective transverse field 
^(12) = hih2/Ji2', the couplings of this new spin to the 
neighbors remain unchanged to leading order. Each such 
cluster c has a moment flc given by the number of initial 
spins in the cluster; when two clusters are combined to 
form a bigger cluster, their moments add: /i(i2) = A'i+M2- 
This procedure is now iterated with the energy cutoff 
n = max{/ij, Jj} of the new effective Hamiltonian being 
gradually lowered. 



In this section, we obtain the long-time asymptotics 
of average imaginary-time autocorrelations in the criti- 
cal region of the RTFIM — we will be using heavily re- 
sults of Ref. ^ refering to sections in that paper by, e.g., 
F Sec. IVB. Our predictions, can be compared with the 
extensive numerical resultal-l available in the literature, 
and this serves as a useful check on the validity of our 
basic approach — the results of s uch a comparison have 
already been summarized in Sec. II G 2. 

We consider the local dynamical susceptibilities 



\{m\af\0)\^Siu-E^)^ 



(61) 



where the sum is over all excited states \m) with excita- 
tion energies Em, and a = x ot a = z. The low- frequency 
behavior of these susceptibilities determines the long- 
time asymptotics of the corresponding imaginary-time 
autocorrelation functions 



C""(t) = «(r)<(0)) 



(62) 



with a = z (local magnetic moment autocorrelation) or 
a = X (local energy autocorrelation); we are consid- 
ering here only the fluctuating (time-dependent) parts 
of autocorrelations and will ignore any constant (time- 
independent) parts (such a constant part in, for example, 
spin autocorrelation represents a non-zero magnetization 
density in the system and is a static property). In the 



21 



following, we simply write Cioc(t) for the local magnetic 
moment autocorrelations C'IJ(r) and Ci^^{t) for the lo- 
cal energy autocorrelations C^j{T) (and similarly for sus- 
ceptibilities). We first obtain (using our basic strategy) 
results for average susceptibilities, which can be conve- 
niently defined as 



[xroc]av(c.) 



(63) 



where L is the size of the system (in the thermodynamic 
limit of i — + CX3 this definition coincides with an ensem- 
ble average over disorder realizations). We also consider a 
semi- infinite chain, > 1, and calculate average dynam- 
ical susceptibilities Xi°' of the boundary spin ci (in this 
case, the average is over disorder realizations). These re- 
sults are then immediately translated to the correspond- 
ing statements about the long-time behavior of average 
autocorrelations. As long as we are only interested in the 
asymptotic behavior of the average dynamical suscepti- 
bilities and autocorrelations, it suffices to use the leading- 
order results for the renormalization of the corresponding 
operators — this is discussed further in Sec. VD. 



We emphasize from the outset that our calculations in 
this section are closely related to the discussion in Ref. ^ 
of static response functions at finite temperature T: Such 
static response properties are calculated by assuming 
that all effective degrees of freedom that are present (in 
the sense of the RG) at energy scale T contribute freely 
to the response at this temperature, while in our calcu- 
lations of the dynamical properties at frequency u> only 
degrees of freedom that are being decimated at scale u) 
contribute to the dynamical response at this frequency; 
even this difference disappears when the dynamical sus- 
ceptibilities are translated to the imaginary time auto- 
correlations, since average autocorrelations at time r ac- 
quire contributions from all frequencies smaller than 1/r. 
The intent here is merely to present a unified approach, 
within the RGoi Ref. ^, to the calculation of such average 
dynamical properties. Also, these calculations, together 
with a detailed physical picture developed in Ref. ^ of the 
phases of the system near the critical point, serve as a 
valuable guide to our intuition in identifying the relevant 
Griffiths regions that dominate a particular response; on 
some occasions in the previous sections (particularly in 
the lAF phase of spin-1/2 chains), such Griffiths argu- 
ments were our only source of information about the be- 
havior of dynamical quantities, and the opportunity to 
compare such suggestive arguments against the results 
of controlled calculations is most welcome. 



1. Average local spin autocorrelation [Cioc]av(T") 

The leading-order renormalizations of the spin op- 
erators are particularly simple: As long as a given spin 
j is active, the operator cr| is renormalized to the "spin" 
operator of the cluster c that the spin j belongs to; 



when this cluster is decimated, the corresponding opera- 
tor renormalizes to zero. 

To calculate [xioc]av('j-'), we run the RG down to energy 
scale fifinai = i-^/S, and rewrite the spectral sum in terms 
of the degrees of freedom of the renormalized problem; 
excitations that contribute to this new sum are clearly 
the (7^ = -l-f excitations of the spin clusters that are 
being frozen into their — —1 states by the transverse 
fields at this scale, and the spectral sum is now easily 
evaluated: 

[Xloc]av(^) = jY,Y.i^c\{m\^mfH^ - Era) 



n(r. 



-Ro{T^)^lfi{T^) , 



(64) 



where we used the fact that all flc spins that are active 
in an effective cluster c contribute identically, and JIqIT) 
is the average magnetic moment of the clusters that are 
being eliminated at scale F. [Note that here, and in the 
following, we simply write instead of T,^/2 to avoid 
clutter in our notation — since we are interested only in 
the leading behavior, the difference is not important for 
our purposes]. 

At criticality (6 = 0), we obtain 



[Xloc]av(w) 



1 



wl InwP 



(65) 



for Lu <ti Qq. For the average spin autocorrelation in 
imaginary time r 3> ^Iq^ we then find 



[Cloc]av(''') 



I In-; 



(66) 



In the disordered phase {S > 0), we obtain 

[xio4.(c.)^<5-^-!i^'^l 



/ziS) 



Inrl 



/z{S) ' 



(67) 



for UJ <^ fls and t ^ fl^ 



Here, we have used scal- 



ing solutions for the off-critical flows to write the an- 
swer for the local susceptibility and have chosen to ex- 
press the power-law in terms of the dynamical exponent 
z{d). From the scaling solution to the RG flow equa- 
tions, we have z~-^ = 2\S\ — this is to be thought of as 
the leading term in a small-(5 expansion for z~^. Written 
in terms of z{6), our result Eq. ( |67| ) is valid more gener- 
ally, and can be understood directly in the simple picture 
of the disordered phase given in F Sec. IVB: The average 
spin autocorrelation at large time t is dominated by the 
(rare) spins that belong to the rare strongly coupled clus- 
ters (Griffiths regions) with low effective "flipping rates" 
(i.e., effective transverse fields) smaller than u; ~ I/r. 
The density of such clusters, which is also the density of 
the most numerous excitations at these low energies, is 
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n(u;) ~ W"'^/^ (this is fixed by the relationship t ^ F be- 
tween length and time scales, that serves as the definition 
of z). Most of these clusters have their effective hipping 
rates between uj and some fraction of lu, and therefore ef- 
fective moments of order | Incjl [since all clusters that are 
being eliminated at a fixed energy scale Q have roughly 
the same magnetic moment proportional to | lnr2|]. Esti- 
mating the contribution of such Griffiths regions clearly 
gives us Eq. ( |67| ) including the factor of | lnr|. 
Finally, in the ordered phase (5 < 0), we obtain 



Int 



^1-2/^(5) 



(68) 



for <C and t ^ ■ In contrast to the case of 
the disordered phase, the interpretation of Eq. (|6|) in 
terms of the picture of the ordered phase presented in 
F Sec. IVA is more subtle. In the ordered phase, the typi- 
cal excitations at low energies w <C fis are classical — they 
are "domain walls" that "break" large clusters apart in 
the places where the clusters are held together by weak 
(effective) bonds of strength of order lu. Such weak ef- 
fective bonds represent the rare large regions (Griffiths 
regions) that are locally in the disordered phase. These 
domain wall excitations are the most numerous excita- 
tions that define the relationship between the energy and 
the length scales and determine the dynamical exponent 
z{6). Such excitations, however, even if they are local- 
ized in the neighborhood of site j, do not contribute to 
Xjji^) since they cannot be "excited" from the (classi- 
cal) ground state by the action of (t| . Excitations which 
do contribute involve much more rare ferromagnetic clus- 
ters that are are flipping back and forth in isolation, with 
flipping rates of order lu ^ 1/r or slower [of course, we 
exclude the macroscopic cluster flipping at a rate of or- 



der 



-cL 



as we are subtracting out the time-independent 



part of the autocorrelation]. In the RG language, these 
are precisely the clusters that are decimated at energy 
scales of order lu, i.e., that happen to have (at these 
scales) an anomalously strong transverse fields of order 
LU [remember that we are in the ordered phase] . A simple 
construction, however, clearly shows that the density of 
such regions is indeed ~ w^/^ , as predicted by the scaling 
solution: For such a cluster to occur we need a ferromag- 
netic segment of length ~ | lna;| (which is not rare in the 
ordered phase) that is isolated (from eventually becom- 
ing a part of the macroscopic cluster) on each side by a 
disordered region of comparable length. Each of the two 
disordered regions is actually a "typical" Griffiths region 
at these energy scales, and the two are required to occur 
much closer together than their typical separation lu^-^^^; 
this explains the appearance of the power 2/z in Eq. (68). 
The factor | lnr| again comes from the typical magnetic 
moments of such ferromagnetic droplets. 



2. Average local energy autocorrelation [Cio(,]av(T) 

We begin by working out the leading-order renormal- 
izations of the operators: When a given spin j is com- 
bined with another spin k into a new cluster (i.e., when 
the strong bond Jjk is being eliminated) the operator crj 
rcnormalizes to {h(^jk)/hj)(7^ji^-^ , where hj is the trans- 
verse field on the spin j before the decimation, ft.(jfe) is 
the effective transverse field on the new cluster (jk), and 
^fjfe) the effective "spin-flip" operator of this cluster 
(this rule ignores a constant term proportional to the 
identity operator, which is unimportant for our purposes 
as we are not interested in the time-independent constant 
piece of the energy autocorrelation function). On the 
other hand, when the spin j is eliminated, the operator 
crj becomes effectively zero to first order in the nearby 
interactions (we again ignore any constants). Iterating 
this, the operator crj is renormalized to {hc/h^^^)a^ if 
the spin j is active in some cluster c with the effective 
held he, and is renormalized to zero if the spin is not 
active; here is the original (bare) transverse field on 
the spin j. 

We now run the RG down to energy scale ilfinai — lu/2 
and rewrite the spectral sum as 



where 



E 

jec 



,(0) 



(69) 



(70) 



In the last equation ^j^^ is over all spins that are ac- 
tive in a given cluster c. Note that at low energies the 
effective cluster field he is only weakly correlated with 
each of the bare flelds ft.!"' ; if the bare field distribution 



is not too broad, we can approximate h 



(0) 



Qq, and 



write gc ~ jic{hc/^o)'^ , where jlc is the moment of the 
cluster c. Doing this clearly misses some nonuniversal 
numerical factor of order one that depends on the bare 
(high-energy) physics. This factor is, in principle, a ran- 
dom quantity that differs from one cluster to another; 
however, this number is expected to be roughly the same 
for all clusters that contribute to the spectral sum at 
low frequencies due to averaging, since such clusters are 
all large, and in some sense similar. Thus, we expect 
that the low-frequency behavior is not affected. [Note 
that we would have been spared this discussion if we 
were to analyze the spectral sum with matrix elements of 
h^P , which is anyway a more natural operator to con- 
sider when thinking of the local energy fluctuations] . The 
excitations that contribute to the spectral sum Eq. ( |69| ) 
correspond to transitions from the cr| = ct^ states to the 
(t| = — (t| states of two (effective) spins j and k that are 
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being combined into one cluster ^(^jk) by a strong bond 
at the energy scale rifinai- Since the log-field distribu- 
tion is broad [we are near the critical point and at low 
energy scales], for such a pair of spins to contribute sig- 
nificantly the transverse field on at least one of the two 
spins involved must be of order lu; thus, we have 

[xroc]av(c^) ^ ^n(r^)Po(r^)i?o(r.)Po(rc.) , (7i) 

"o 

from which we immediately read-off our results: 
At criticality, for w <C fio a-nd r 3> , we obtain 



[Cfoclav (r) 



|lnw|4-<^ ' 

1 

r2|lnT|4-<#' 



(72) 



Away from the critical point, both in the disordered 
phase and in the ordered phase, we obtain 



[Cfoclav M ~ I'^l^ ''' j+2/}{S) ' 



(73) 



for ui <C and r 3> ^J^', in the last formula we again 
used z{6) as a more physical parameter characterizing the 
Griffiths phase at a given 5. The off-critical energy au- 
tocorrelation function thus behaves similarly in the two 
phases, as expected from duality. It is again possible to 
interpret these results in terms of the statistical proper- 
ties of appropriate rare regions that dominate the aver- 
age energy autocorrelation at long times. As the results 
(and their interpretation) are identical in either phase, we 
sketch only the interpretation on the ordered side: As we 
have already noted, for a region to have significant energy 
fluctuations at the frequency scale of order w, it must con- 
tain two adjoining segments both having a characteristic 
energy of order to — a predominantly disordered segment 
(in the RG language, an effective bond with J ^ uj across 
this segment) and a predominantly ordered segment (in 
the RG language, a cluster with h ^ cu). Clearly, the 
predominantly ordered segment with effective transverse 
fleld ~ uj can exist only if it is also isolated on the other 
side from the rest of the system by another predomi- 
nantly disordered segment having the same characteristic 
energy scale. This situation has already been analyzed 
in the context of the spin autocorrelations in the ordered 
phase, and clearly one recovers precisely Eq. ( [73| ) from 
such an analysis. 



3. Autocorrelations of the boundary spin 

So far, we have calculated average autocorrelations for 
the spins in the bulk; calculations for the first spin cri in 
a semi-infinite chain j > 1 proceed analogously, and we 
will simply state the results. 



For the spin autocorrelation, we find 

^(o|r^) Po(r^)flo(r^) 

[XlJav(t^) = 



(74) 



here E{l3\T)dp is the probability that ai survives to scale 
r and is in a cluster with \n{^/h) = /?; such properties 
of the boundary spin are fully characterized in F Sec. V. 
At criticality we obtain 



[Xllav i^) 

[Cilav (r) 



1 



I In w p 
1 



llnri' 

while away from the critical point, we get 



(75) 



[Xllav (^) 



(76) 



As in the bulk case, we can interpret the average off- 
critical spin autocorrelation Eq. (^6|) in terms of the rare 
instances that dominate this average. In this case, the 
corresponding rare regions must start at ai — this ex- 
plains the absence of a | lnr| factor in Eq. (|7^) compared 
to the bulk resuhs Eqs. 1^) and (||). The only other 
difference is that in the ordered phase we do not need 
to isolate the ferromagnetic droplet (containing ai ) from 
the left. 

For the energy autocorrelation we find 



[x?]av(w)^^Po'(rc.)i?o(rj 



(77) 



At criticality we obtain 



[CllavW 



lntj|3 
1 



'^^^ ' r2|lnr|3 ' 
while in the disordered phase 

[Xll^ (c) ^ <53^i+2/.(.) 

and in the ordered phase 

[x!L(c.)^|<5|3c.i+™, 

|<5|3 



(78) 



(79) 



(80) 



Thus, the average off-critical energy autocorrelation of 
the boundary spin differs from that of the bulk spins in 
exactly the same way (and for the same reasons) as in 
the case of the spin autocorrelation. 
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C. Dynamic structure factor of the spins 

Let us now briefly consider the dynamic structure fac- 
tor S^^{k,uj) defined as 

5- ik,u;)^jJ2\ (m| J2 ^'"^^ '^j m'Siu;-E^). (81) 



S^^{k,uj) characterizes the spatial structure of the exci- 
tations at energy uj. Proceeding as before, we find 

S^^ik, u) ^ !^R,(T^)\iI^{r^) , (82) 



where \^o{k)\^{T) is the average modulus squared of the 
effective magnetic moment at wavevector k for the clus- 
ters that are being eliminated at scale F; for a given 
cluster c, this effective moment is defined as /Ltc(fc) = 
X^jGc ■ dynamic structure factor can also be 

written in terms of the function D{P, x\Tt^) defined in 
F Sec. IIIB4; we have 



Dio,k\r^) 



OJ 



(83) 



where D{0, k\r^) is the Fourier transform of D{0,x\r^) 
at wavevector k. We have not attempted to analyze 
D(P,x\Ti^), even though a detailed characterization is 
likely to be possible (see F Sec. IIIB4). Instead, we will 
only analyze the behavior of the dynamic structure factor 
in some limiting cases using the scaling picture. 

First, consider the system at criticality. Fix wavevector 
k. Then, for F ^ F^ = 1/^/k the effective cluster mo- 
ments at wavevector k "add coherently" (more precisely, 
the real parts of the effective moments of the clusters 
that are being combined into bigger clusters are of the 
same sign) and therefore scale as F*^. At scales F » F^, 
the effective moments at k "add incoherently" (the real 
parts of the moments being combined can be of any rel- 
ative sign) and therefore scale as F^(F/Ffe)'^='y'", where 
0sym = (1 + \/5)/4 is the growth exponent for the cluster 
moments distributed symmetrically around zero (see Ap- 
pendix of Ref. H). Thus, we arrive at the following scaling 
form for the dynamic structure factor at criticality 



S^'{k,uj) 



p20 

^ UJ 



^{kTl 



(84) 



where $(a;) ~ const for a; ^ 1 and <i>(x) ~ for 
a; 1. We cannot, however, address the regime /cF^ ^ 1 
by such a scaling analysis. 

Now, consider the system that is not critical, either 
in the disordered or in the ordered phase, in the regime 
F^j ^ Ts- The length and the magnetic moment of a 
cluster that is eliminated at scale F ^ Ts are sharply 
defined: lo{T) = q(F/F5)F2 + Oiy^TfTsTj) and Mo(fc = 
0,F) = c^(F/F5)F^ + Oiy^TJTsVf), where q and are 
numerical constants of order one. Such a cluster has some 



internal structure on the length scales below the correla- 
tion length 6~'^, but "looks" fairly uniform on larger 
length scales. Then, for the wavevectors k 5^ we have 



^4-20 
fc2 



l + cos(Q/cF/|(5|)e-'='='r/|5|^ 



(85) 



note the oscillatory behavior at k ~ |<^|/F due to the 
"sharpness" of the lengths Iq of the clusters; the gradual 
suppression of this oscillatory behavior at larger wavevec- 
tors comes from the uncertainty in Iq, which is much 
smaller than Iq itself. To obtain the dynamic structure 
factor S^^{k,uj) we simply need to multiply this "cluster 
structure factor" by the density p{ijj) of such clusters at 
energy uj: p{u)) ~ J3^^i-i/z((5) ^^le disordered phase 
and p{uj) ^ J'^/w^^^/^^'^^ in the ordered phase. 



D. On the validity of the results 



Our analysis in Sec. Ill E of the effects on the (calcu- 
lated) dynamic spin structure factor of higher-order cor- 
rections to the spin operator renormalization rules in the 
spin-1/2 antiferromagnetic chains can be carried over to 
the quantum Ising spin chain as well. Such higher-order 
corrections will only renormalize the effective magnetic 
moments of the remaining clusters by numerical factors of 
order one — this is also discussed in detail in F Sec. VIA. 
We do not repeat such analysis here. Instead, we consider 
the effect of such corrections on the average autocorrela- 
tions of the boundary spin ai. This is somewhat more 
tractable than the bulk case, and we can actually include 
these effects in an explicit calculation. Our results are 
not surprising: these next-order corrections only affect 
the values of some non-universal prefactors, and have no 
effect on the time dependence of the autocorrelations in 
the long-time limit. 

The first-order renormalization of the boundary spin 
operator af is discussed in detail in F Sec. VB in the 
study of the end-point magnetization, and also in Ref. ^ 
in the study of the end-to-end correlations. Following 
these references, we write the component of the boundary 
spin on the left-most remaining cluster ai as (crfaf) — 
e^^. Then, for the special family of scaling solutions, the 
distribution of the parameter A at scale F is described by 
the probability density 



C{A\r) = ^^SiA) 
>o(,i I) 



Fo(F,)-Po(F) 
^^o(r/) 



Po(r)e- 



-Po(r)A 



(86) 



where the first term is simply the probability that ai 
survives down to the scale F from the initial scale F/ (we 
already considered this probability in our zeroth-order 
calculation). The average autocorrelation is then given 
as 
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[Xl]av(w) 



LoPoiTi) 



i^o(r/)-Po(r^) 



Po(r^) 



(87) 



where we have kept aU the numerical factors within the 
RG; the first term comes from the instances when the 
boundary spin survives down to energy scale and is 
precisely our zeroth-order result Eq. ([74|). In all cases of 
interest Po{Tu) <C 1, and the second term is thus roughly 
the first term times Po(rj). Now, l/Po{Tj) is the width 
of the initial distribution of the logarithms of interac- 
tions, and is therefore some number of order one — the 
two contributions are thus of the same order. However, 
the wider this initial distribution is, the smaller the sec- 
ond term is relative to the first. The origin of the second 
term is easily understood: the boundary spin tJi may 
be decimated with some finite probability at an early 
stage of the RG; in such a case cti will have a signif- 
icant component on the nearest surviving spin, which 
may now be considered as a boundary spin in the new 
problem with a somewhat smaller energy cutoff; the av- 
erage large-time (low-frequency) autocorrelation of this 
new boundary spin will be at least of order Eq. (jTJ) from 
the instances when this new spin survives to the energy 
scale uj. The above calculation explicitly shows that the 
average dynamical response of the spin ai is dominated 
entirely by the instances when one of the spins within 
some small (of order one) distance from the boundary 
survives down to low energies of order oj. 

Thus, we see that our general arguments for the ef- 
fect of higher-order terms in the operator renormaliza- 
tion rules are borne out by this detailed calculation for 
the boundary spin autocorrelation. This suggests that 
our approach gives asymptotically exact results (apart 
from non-universal prefactors) in the bulk case as well. 



VI. A DISCUSSION OF T / PROPERTIES 

So far we have calculated various dynamical and trans- 
port quantities at T = 0. These results are clearly valid 
even at T 7^ so long as the probe frequency u sat- 
isfies T UJ. Unfortunately, it is not straightforward 
to generalize our calculations to the complementary low- 
frequency regime (w <C T) dominated by thermal effects. 
There is, however, one exception. As mentioned earlier, 
the spin- 1/2 XX chain is equivalent to a model of spinless 
fermions with random nearest neighbor hopping and zero 
chemical potential. It should come as no surprise that the 
free-fermion nature of this problem allows us to straight- 
forwardly calculate some dynamical and transport prop- 
erties at small non-zero temperature. We begin this sec- 
tion by formulating a fermion analog of the RG procedure 
used for the spin chains. We will then use this RG ap- 
proach to work out the low-frequency, low-temperature 
dynamical conductivity and the zz component of the dy- 
namic structure factor for the spin-1/2 XX chain without 



any restriction on the relative magnitudes of u and T (the 
calculation of the perpendicular component of the struc- 
ture factor for cj < T is much more complicated, and 
we will only be able to discuss its qualitative behavior). 
Naturally, these results are not at all generic, relying as 
they do on the free-fermion nature of the problem. On 
the other hand, a weak coupling, which corresponds 
to the nearest-neighbor interaction between the fermions, 
is strongly irrelevant in the RG sense at the free-fermion 
point, and the system flows to the non- interacting point. 
Since this non-interacting limit is singular as far as finite 
temperature transport properties are concerned, we have 
here an example of a "dangerously irrelevant operator" , 
and the important physical question is how this weak, ir- 
relevant interaction affiects the T 7^ transport near the 
non-interacting point .EJ This is what we turn to at the 
end of this section. 



A. Free fermion RG 



The free-fermion problem Ti. — J^j 



^3+ 



has been the subject of extensive investigation in the past 
using a variety of techniques (see, e.g., Ref. ^ and ref- 
erences therein). For our purposes, it is most convenient 
to introduce a RG procedure analogous to the singlet 
RG used in the spin problem. We formulate this proce- 
dure directly in terms of the corresponding single-particle 
Shrodinger problem H = V . t,{\j){j + 1| + |j + l)(j|); 
this RG is, for the case of the Hamiltonian above, es- 
sentially just an efficient way of (approximately) diago- 
nalizing random symmetric tridiagonal matrices with ze- 
roes on the diagonal. We begin with the observation that 
the particle-hole symmetry of the problem causes eigen- 
states to occur in pairs, with energies ±e. The strong- 
randomness RG proceeds by eliminating, at each step, 
such a pair of states with energies at the top and bottom 
of the band: One finds the largest (in absolute value) 
hopping amplitude in the system, say ^2 connecting sites 
2 and 3; this defines the bandwidth Hq = 2 x max{\tj\} 
of the original problem. If the distribution of the tj is 
broad, the symmetric and antisymmetric wavefunctions 
living on these two sites will be good approximations to 
eigenstates with energies ±r2o/2, as ^1/3 will typically 
be much smaller in magnitude than ^2- The couplings 
can then be treated perturbatively, and eliminat- 
ing the high-energy states living on the sites 2 and 3 
results in an effective hopping amplitude ti = —tit'i/t2 
between the neighboring sites 1 and 4. More precisely, 
in the effective Hamiltonian that describes the remain- 
ing L — 2 states, the block 1-2-3-4 is represented as 
7Yi„4 = fi(|i)(4|-h|4)(i|), where the states |i) and |4) are 
essentially the original |1) and |4) states up to 0(^1/3/^2) 
corrections. This rule is essentially identical to the rule 
for the singlet RG at the XX point, as the additional 
minus sign can be 'gauged away' in the nearest-neighbor 
model in one dimension; we will, in fact, only keep track 
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of the absolute values of the tj. The distribution of \tj\ 
in the renormalized problem with bandwidth Q will thus 
be the same as the renormalized distribution of J-^ at 
cut-off scale ft in the singlet RG for the spin problem. 
The analysis of the asymptotic validity of this approach 
thus carries over unchanged from the singlet RG. 

This procedure can therefore be iterated to reach lower 
and lower energies; at each stage we trade in our current 
problem for a new problem defined on two fewer sites. 
This new problem will have the same low-energy eigen- 
values as our original problem. However, evaluating ma- 
trix elements of operators between two low-energy states 
requires some care, as the states |j) in terms of which 
the renormalized problem is written are different from 
the states \j) of the original problem. As in the singlet 
RG, this is best handled by renormalizing the operators 
as we go along, so that the matrix elements of the renor- 
malized operators between the states of the new problem 
are the same as the matrix elements of the bare operators 
between the corresponding states of the original problem. 
This allows us to calculate various dynamical properties 
by evaluating the corresponding spectral sums exactly as 
in the spin language. At T = 0, this amounts to noth- 
ing more than a restatement in terms of the fermions of 
our previous calculations. The new language, however, 
has one important advantage: thermal effects are easily 
incorporated into this framework, essentially because the 
non-interacting nature of the problem is made explicit. 
[We emphasize again that the RG finds all cigenstates of 
the free-fermion problem. The corresponding statement 
can also be made in terms of the singlet RG on the XX 
spin chain: when eliminating a pair of spins 2 and 3 the 
effective Hamiltonians in all sectors (corresponding to the 
states |so), \to), and \t±i), of the pair) are identical up 
to a sign of J^^ in the |t±i) sector.] Moreover, it is now 
possible to address questions specific to the free-fermion 
problem that do not have a natural analog in the spin 
language (see Appendix P) . 

Finally, we note in passing that this RG procedure 
can be generalized to analyze other particle-hole sym- 
metric free-fermion problems in one and two dimensions 
(which are not immediately equivalent to any quantum 
spin problem) as well as analyze the general properties of 
the Bogoliubov-de Gennes equation for quasiparticles in 
a one dimensional superconducting wire in the absence of 
spin-rotation symmetry (the results of such an analysis 
will be published separately). 



B. T 7^ dynamics and transport at the XX point 

Let us begin by working out the full T and lu depen- 
dence of the dynamical conductivity, Eq. (^ , at the free 
fermion point. Our first task is to work out the rules that 
govern the renormalization of the current operators T(j). 
Assume once again that the hopping element t2 has max- 
imum magnitude. We wish to work out what operators 



we should use in place of T(l), T(2) and T(3) when we 
renormalize down to lower energies by eliminating the 
corresponding two states at the top and bottom of the 
band (the other current operators to the left and right 
of this segment will be unchanged to leading order by 
this elimination). An explicit perturbative calculation 
immediately yields f(2) = t(l/3) = i?i(|i)(4| - |4)(i|); 
this is completely analogous to the rule obtained in the 
spin representation, and as before, we will call this oper- 
ator T(l) for consistency of notation. 

As we carry out the RG and reduce the bandwidth by 
removing states from the top and bottom of the band, 
the above result implies that T(j) renormalizes to 

^^ZjT(j), where j now labels the sites of the renor- 
malized problem with bandwidth Q, and the Ij are the 
lengths of the renormalized bonds in this problem. With 
this in hand, we run the RG until the bandwidth is re- 
duced to Ofinai = uj and rewrite the spectral sum Eq. ( [55| ) 
as 

^1,^2 j 

x[f{e^,)-f{e^,mcj-e^,+e^,). (88) 

Because of the extremely broad distribution of the ij , the 
dominant contribution to the sum Eq. ( |8^ ) comes from 
transitions between the two members, one at the bottom 
and the other at the top of the renormalized band, of 
each pair of states that is being eliminated at this en- 
ergy scale. The matrix element for this transition is just 
Z|a;|/2, where I is the length of the hop in question. In 
the thermodynamic limit, we thus have 



a'icj) = [fi-cj/2) - f{cj/2)] X 

niTui) 2;2 



dl dCio'^P P{C, 1\V^)5{lo - Loe-^) 



sinh(c./2r) ^, ^ 
2cosh2(tj/4T) ^-^^ ' 



(89) 
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FIG. 7. A plot of the frequency dependence of o-'(cij) at low 
T 7^ at the XX point. Note that this result is expected to 
break down below a frequency scale l/rtr <^T w hen a small 
but non-zero interaction is turned on (see Sec. VIC). 
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which is the leading behavior for w, T <C rJo- This result 
smoothly interpolates between the logarithmic frequency 
dependence seen earlier for 3> T and the limiting form 
a'{uj) ~ ujln{no/uj)/T valid for uj <C T — a plot of this 
frequency dependence is shown in Fig. 0. 

Let us now turn to the T ^ spin dynamic structure 
factor at low frequencies in the vicinity of k = ir/a. In 
the single-particle language, the spectral representation 
for S^^(k, uj) reads 

x[f{e^,)~fie,,mu;-e^,+e^J. (90) 

Here S^{k) = ^^ 8^(^)6'*^^^ is the Fourier transform 
of the position-dependent matrix operator S^(j); in the 
real-space basis S^(j) = n{j) — 1/2, where n(j) ~ 
This spectral sum can also be evaluated within our RG 
approach. The leading-order operator renormalizations 
in this case are, in complete analogy with the spin prob- 
lem, very simple: each S^{j) remains unchanged unless 
a state living on site j is eliminated, in which case S^{j) 
renormalizes to a multiple of the identity. As before, we 
run the RG till the bandwidth is reduced to fifinai = 
and do the spectral sum with the renormalized operators 
in the new problem. This renormalized sum may be eval- 
uated by again recognizing that it is dominated by tran- 
sitions between pairs of states with energies ztuj/2 that 
live on pairs of sites connected by 'strong' hopping ampli- 
tudes (of magnitude uj/2) in the renormalized problem. 
The corresponding matrix element is just (1 — e'*^')/2, 
where I is the length of the hop in question. Counting 
the contributions exactly as in our zero-temperature cal- 
culations, we thus get 

S''ik,Lu) = X S^UiKi^) ■ (91) 

Thus, we see that S^^ is essentially unaffected by ther- 
mal fluctuations at the XX point; in particular, the low- 
frequency divergence is not cut-off by temperature effects 
even when lu <^ T. 

A similar analysis can clearly be performed in the XX- 
RD phase. Again, both a' (to) and S^^{k,u!) at T > are 
simply given by the corresponding expressions at T = 
multiplied by simple functions of io/T, exactly as in 
Eqs. ( |89| ) and (|9|). Thus, though temperature effects 
are simple to work out at the XX point and in the XX- 
RD phase, the results are rather special due to the free- 
fermion character of the problem. 

C. Going beyond the free-fermion results 

What happens when we turn on the nearest-neighbor 
interaction? This is the question we need to address next. 

Let us first consider the effects of small couplings 
added on to the XX model. The analysis of Ref. ^ shows 



that this term is irrelevant in the RG sense; the typi- 
cal value of J^/J-"- at log-cutoff F scales as (J^/J-*-) ~ 
uo exp(— cF-^), where c is an 0(1) constant, (p is the 
golden mean (1 + %/5)/2, and uq is the typical value 
of J^/J-^ in the microscopic model. A useful way of 
thinking about the low-frequency behavior of the con- 
ductivity is as follows: Imagine running the RG till the 
cutoff n ^ T. In this renormalized problem the typ- 
ical {.P/J^) uoexp(-cF^), where Ft = ln(r2o/r). 
In the fermion language, this is the typical value of the 
ratio of the nearest neighbor interactions to the hop- 
ping amplitudes. In this renormalized problem, a naive 
Fermi's Golden Rule estimate of the corresponding in- 
elastic collision rate due to interactions gives I/tcoIi ~ 
Mgrexp(— 2cF^). This gives us a frequency scale below 
which our T 7^ free-fermion results are expected to 
break down as a result of the residual interaction effects. 

Unfortunately, we are unable to do a controlled cal- 
culation that determines the transport properties in the 
frequency regime uj < 1/tcou. The best we can do is 
to work out what a naive scaling argument would pre- 
dict for the d.c. limit of the conductivity. The basic idea 
is as follows: The collision rate may be converted into 
a corresponding dephasing length Lcou by appealing to 
the activated scaling that is a characteristic of our prob- 
lem. This gives Lcoii ^ lii^(''coii) ^ r^^- This is the 
length scale beyond which quantum coherence is lost due 
to inelastic collisions. Now, we can imagine breaking 
up the system into blocks of length Lcoii- A d.c. cur- 
rent / passing through the system will see a chain of 
resistors corresponding to these blocks — the resistance 
values of each of these blocks is simply given by the 
T = Landauer resistance of the corresponding system 
of length Lcoii. The voltage developed across a system 
of total length L will therefore he V — I J^fLi''"" — 
ILRa,y{Lco\\)/Lco\\- Since -Rav(-Z^coii) ^ e''^'^'"'" (where ci 
is an 0(1) scale factor), the d.c. conductivity works out 

to be cTd.c. ^ Leoue-^i^-" ~ ln2^(r!o/T)e-^''"'"("''/^). 
Note that in the absence of interactions, we had earlier 
found cr'{Ld) as w at T > — our scaling ar- 
gument implies that interactions render this conclusion 
invalid. Unfortunately, while this scaling argument is cer- 
tainly plausible, the question of the true low-frequency 
limit can only be settled by a controlled calculation in 
the regime tu <§; 1/tco1i, which is beyond our current ca- 
pabilities. 

The above arguments also suggest that S^'^{k,oj) will 
deviate from the T ^ free-fermion result for w < 1 /tcoU. 
In particular, one expects that the w — > divergence of 
S^^ would be cut-off below this frequency scale. Sim- 
ilar behavior is also expected of S'"' , but again, what 
is really needed is a controlled calculation as opposed to 
a scaling argument. Note also that we expect something 
different at the XXX and XXZC quantum critical points: 
since the theory at these critical points already includes 
interactions, one expects that I/tcoh ~ T, and the relax- 
ational behavior characteristic of an interacting system 
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at finite temperature will set in for w ~ T, in contrast to 
the behavior in the vicinity of the XX point. 



VII. PROSPECTS FOR EXPERIMENTAL TESTS 

Previous experimental work on one-dimensional 
random-exchange Heisenberg antiferromagnetic spin 
chains has characterized the dynamics of these systeim 
in terms of the observed NMR l/Ti relaxation rate,E3 
and ESR relaxation rates and linewidths.Eil 

As far as the NMR measurements are concerned, our 
calculations are unfortunately not directly relevant to the 
experimental measurements of 1 /Ti. This may be seen as 
follows: In the usual case of a pure, translationally invari- 
ant system, 1/Ti is directly related, by Fermi's Golden 
Rule, to the local dynamic structure factor ^loc evaluated 
at frequency w equal to the nuclear resonance frequency 
"/nH , where is the nuclear magnetic moment and H 
is the external field. In a random system, with a broad 
variation in the value of S'ioc('^), the following question 
immediately arises: what measure of the distribution of 
Sioc{oj) does the experimentally measured l/Ti{H) re- 
flect? 

Now, we have seen that the average 5*100(1^) diverges 
strongly as w — > at T = 0. Naively, one might have 
thought that this would imply a corresponding diver- 
gence in 1/Ti at small H, at least when T <^ H. How- 
ever, the divergence in 5ioc(w) comes from a few very rare 
sites which give a very large contribution. Clearly, the 
observed l/Ti will be completely insensitive to this effect, 
since all that will happen is that a tiny fraction of nu- 
clear spins (in the neighborhoods of those rare electron 
spins that have significant spin fluctuations at the fre- 
quency to = "/nH) will flip almost instantaneously, while 
the rest of the nuclear spins will have an extremely small 
probability to flip, and this is what will be reflected in 
the spin relaxation experiments. In this sense, it is the 
typical value of Sioduj) that is more relevant for com- 
parisons with NMR 1/Ti data. A typical nuclear spin 
will in fact have essentially no spin fluctuations to cou- 
ple to at frequency u = ^nH — it can therefore relax 
only by paying an activation energy that is set by jeH 
(where 7e 3> is the electron magnetic moment) since 
the external field acts to freeze out all modes below this 
energy scale in most of the system (with the exception 
of the rare regions alluded to above). The experiments 
actually do see activated behavior for 1 /Ti at finite tem- 
perature. However, the activation gap seems to scale as 
A ~ H^-^ — the rough argument above of course cannot 
explain this non-trivial H dependence of the observed 
activation energy. 

Our second remark relates to th£ ESR linewidth mea- 
surements of Clark and Tippie.Ell Here, again, our re- 
sults do not address the experimentally relevant ques- 
tions. This is because all our calculations for the XXX 
case are done within the context of the simple Heisenberg 



exchange Hamiltonian, while the observed linewidth in 
the experiments is determined by other effects such as 
dipolar interactions or anisotropy. 

Inelastic neutron scattering experiments, on the other 
hand, if they can be done on these systems, provide a di- 
rect testing ground for our predictions. We conclude with 
some remarks on the relevance of our calculations of the 
dynamic structure factor to such experiments. First of 
all, note that we considered randomness in the exchanges 
only, with the spins themselves assumed positioned on 
regular lattice sites; thus, our results are restricted to 
compounds with chemical disorder in exchanges. It is 
clear that small randomness in the positions of spins (e.g., 
due to thermal fluctuations) will result only in some fur- 
ther suppression (by the standard Debye- Waller factor 
at wavevector k) of the features relative to an overall 
background. In the dimer phase, a possible difference 
in the lengths of even and odd bonds will result only in 
some phase factor in the cosine of Eq. (28). Also note 
that the non-magnetic neutron scattering from such spin 
chains will actually be suppressed near k = tt/o, and this 
may facilitate a possible experimental observation of the 
predicted features. We caution, however, that while it 
would be extremely interesting to see the sharp oscilla- 
tory structure predicted in the Griffiths phases, this may 
be difficult to achieve without going to very low tem- 
peratures and energy transfers. Regarding transport, we 
hope that our results will motivate experiments to probe 
the spin conductivity in these systems. 
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APPENDIX A: FINITE-SIZE SCALING 
FUNCTION FOR THE CONDUCTIVITY 

Consider a finite chain with an even number of sites 
Nj = L + 1, where L is the length of the chain, and with 
free boundary conditions (a similar analysis can be car- 
ried out for chains with an odd number of sites, and also 
for chains with periodic boundary conditions). We want 
to calculate the real part of the dynamical conductivity 
averaged over the distribution of bond strengths in the 
limit of low frequencies and large L. We work this out 
for the XX chain; the result in the presence of cou- 
plings will differ only in the values of some non-universal 
scale factors so long as the system does not develop Ising 
antiferromagnetic order in the thermodynamic limit. To 
proceed, we need to keep track of the joint distribution 
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at scale F of the number of remaining spins N, the A'' — 1 
couphngs Q, and the corresponding bond lengths U. In 
a finite system, the couplings become correlated due to 
the constraint imposed by the finite length pf the sys- 
tem. However, following Fisher and Young,L3 we note 
that the couplings remain 'quasi-independent', and can 
be described in terms of the infinite-chain distribution 
P{C, ^|F) exactly as in Ref. |3^. More precisely, if we also 
keep track of the lengths Ip and Ir of the 'dead' regions 
(consisting of singlet pairs formed at higher energy scales) 
at the left and right ends of the chain, then a distribution 
of the form 

dProb[7V; Ci, ^1 • • • CN-ijN~i;lF,h}L, F] = 

= aN{L\T)PiCiJi)dCi . . . P{CN^i,lN~i)dCN~i 
X C{If)C{Ir)5i^+...+i^_^+Ip+i^^l (A1) 

for even N > 2 has its from preserved under renormal- 
ization if aAr(L|F) is independent of N with 

^^ = 2P„{T) = 2jdlP{0,l). (A2) 

Here, P{C, Z|F) satisfies the same flow equation as in the 
infinite chain, and C{l\r) satisfies 



dC 



/ p{c,-)dC-PoC. 

Jo 



(A3) 



In the above, the F dependence is left implicit, and 
/(■) *i is used to denote a (discrete) convolution in 
the length variables. For clarity, we work explicitly with 
discrete lengths, with Ip and Ir even, and li odd integers; 
this is clearly preserved under the RG. 

We start the RG with JIq = l, T/ = 0, the initial bond 
distribution P(C|F/) = (this corresponds simply to 
choosing the initial to be uniformly distributed in the 
interval [0,1]), k — 1, Ip — Ir — 0, and Nj = L + 1; with 
initial distributions P{(,l\Tj) and C{l\r) normalized to 
unity, the normalization factor is a(F/) = 1. The dynam- 
ical conductivity is now given by 



(A4) 



with a(F) = (F + 1)^ (for our specific choice of initial 
conditions), and 



L+l 

A(i|F) = (7V-1) mii)ii 

N = 2 h,l2,...lN-ldF,lR 



P{C2,h)dC2. 



PiCN-l,lN-l)dCN-l 



(A5) 



where the sum is over even N. 

Now, multiplying A(i|F) by e~^^ and summing over 
odd L > 1 — i.e., doing a (discrete) Laplace transform in 
L — removes the constraint on the lengths, and we find 



Aiy\T)=C'iy)Q{y) 



(1-T2(y))2 



(A6) 



where Q{y) and T{y) are respectively the Laplace trans- 
forms of P{0,l)P and f^P{C,l)dC- Thus, we can 
straightforwardly work out My), given P(C, y) and C{y). 
Using the results of Refs. ^ and |[ we can write the 
following explicit expressions for these two functions: 



P(C,2/|F) = r(y|F)e-«<«in, 
uiO\T)uiy\Ti) 
^(^l^) = .(,|F).(0|F,)^^^l^^)' 



(A7) 
(A8) 



where u{y\T) = D{y) coth[D{y){r + C{y))] and r(y|F) 
D{y)/ smh[D{y){T + C{y))]. The functions D{y) and 
C{y) depend on the initial distribution P{C, y|F/), and in 
our case are given by D{y) 
ln(l 



VI - e-2y and D{y)C{y) = 
Also, C{y\rj) = 1. 



y + ln(i + Vl - e-2w) 

With this in hand, it is a relatively simple matter 
to work out A(L|F), L odd, by performing the inverse 
Laplace transform: 



-1 pc+iir/2 

A(L)^- A{y)ey'^dy. 

J c-iiT/2 



(A9) 



In Sec. HID 5, we evaluated this integral numerically to 



compare the RG predictions with the results of the exact- 
diagonalization studies (note that in the main body of the 
paper we didn't make a distinction between Ni and L, 
since it is irrelevant in the thermodynamic limit; in the 
more detailed notatio n of this Appendix, our numerical 
results of Sec. HI D 5 are for system sizes Nj = 128 and 
256). ' 

In the scaling limit F 1, L ^ 1, the integral Eq. ( [A9| ) 
for ^(L|F) is dominated by small y and can be approxi- 
mated by 



A{L\T) ^ 2LT-' A{y\T) . 



(AlO) 



where LT~^ denotes the inverse of the continuous Laplace 
transform. Moreover, in this limit, A(y) may be worked 
out using the following scaling forms for C{y) and P{C, y)- 



F^/2^coth[F^/2^ 



(All) 



P(C, y\T) = . e^Cv^cothFv^] (^12) 

smh[I y/^y\ 

Putting everything together, we can now write A(L, F) = 
Tf{T'^/L), which immediately implies a scaling form for 
the conductivity: cr'(w,L) = T^Q(r'^/L). Thus, we 
see that the dynamical conductivity in a finite system 
satisfies a scaling form that reflects the activated dy- 
namical scaling at the XX fixed point. Note that while 
the scaling form holds more generally, the values of the 
non-universal scale factors that we have used are spe- 
cific to our choice of initial distribution. Analyzing the 
behavior starting with an arbitrary initial distribution 
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allows one to relate these non-universal scale factors to 
the properties of the initial distribution under the as- 
sumption that 'bad decimations' early in the RG do not 
affect these values. Such an analysis allows us to write 
a'{uj,L) = ly\Yi{VlQ/uj)Q (l^hi' {VIq/lo) / L) , with the mi- 
croscopic energy scale ^Iq and the microscopic length 
scale Iv precisely as defined in the main text. Moreover, 
it is clear that the same scaling function also describes 
the low-frequency dynamical conductivity in a large but 
finite system even in the presence of interactions as 
long as the system is in a Random Singlet state — only 
the values of the non-universal scale factors are expected 
to change. 

While it is possible to calculate the full scaling function 
Q{x) by a detailed analysis of the inverse Laplace trans- 
form, we will confine ourselves here to working out Q{x) 
in two limiting cases: For x <^ 1, 8(x) — 7/180 (which 
correctly reproduces the infinite-size result, as it must), 
while in the limit a; ^ 1 we have Q{x) = e~^/^/\/27rx. 
This is t he result used in the Griffiths argument of 
Sec. fnD. 



APPENDIX B: LOW-ENERGY 
EIGENFUNCTIONS OF THE 
ONE-DIMENSIONAL RANDOM HOPPING 
HAMILTONIAN 

We have already noted on several occasions that sub- 
leading terms in our evaluation of low-energy properties 
via the strong-randomness RG can be treated consis- 
tently within the RG picture, and are naturally related 
to the ground state correlations in the system. Here we 
discuss more explicitly one example of this connection, 
for the case of the XX model; our discussion is couched 
completely in the free-fermion la nguag e and uses the 



fermionic RG introduced in Sect. VIA — this will allow 
us to work in a more familiar context and to connect 
with the known results for the fermionic problem. 

Our example relates to the structure of the low- 
energy eigenstates of the one-dimensional random hop- 
ping ^ Hamihonian H = Y.^'^ ( | J + 1) (j | + | j) (j + 1 1 ) ; 
for simplicity, we work with free boundary conditions. 
We note at the very outset that the average Greens 
functions [G(a;, a;'; £')]av have been calculated exactly in 
Ref. using SUSY techniques. Moreover, a great deal 
is known about the statistics of the low-energjLJvavefunc- 
tions, via a sophisticated real space analysis. Here, we 
present a simple (if somewhat rough) picture of the wave- 
functions which allows one to easily reproduce the previ- 
ously published SUSY results, as well as provides us with 
the following rough picture of a typical pair of wavefunc- 
tions at energies ±e: such a pair is typically "localized" 
on two sites i and j separated by a distance of order , 
= U{\i)±\j))-^. . . with TV of order 1, with the wave- 
function amplitudes decaying as e^"^^ away from the two 
sites for r <• F? and as e~^ for r > F?. Moreover, as 



we shall see below, the above picture is intimately related 
to corresponding structure in the zero-energy wavefunc- 
tion. 

We start by considering this zero-energy wavcfunction. 
If the number of sites L is odd, the Hamiltonian H 
has a (unique) zero-energy eigenstate '0, and this cigcn- 
state can be written down immediately as '!/'2n = and 
"0271-1-1 = H"^j(-i2j_i/t2i)'0i- Thus, the logarithm of 
the wavefunction In \^2n+i \ performs a random walk (see 
Fig. H). While the relative amplitudes of the wavefunc- 
tion are thus understood directly, we are interested in 
the properties of the normalized zero-energy wavefunc- 
tion 5", that is, we also need to understand normalization 
of the eigenvector ip. Such properties of the normalized, 
wavefunction ^E" were analyzed by Balents and FisherCJ 
by reducing the calculation of correlators to the Liouville 
quantum mechanics. Here we suggest a more direct ap- 
proach by noting that normalization of the -0 for a given 
realization can be achieved by simply setting the largest 
amplitude V'max — max{'02n-i-i} to a number of order one. 
More precisely, if we set V'max = 1 , then for any g > the 
norm J\f2q = IV'iP'^ is a random quantity of order 

one with a well-behaved ( "narrow" ) distribution of width 
of order one. This can be seen directly by, for example, 
calculating first moments of the distribution of Mq us- 
ing the properties of such a random walk near its global 
maximum (absorbing boundary). With this in hand, we 
can characterize the zero-energy wavefunction essentially 
quantitatively even without the knowledge of the right 
normalization constant, which is always of order one. 

Let us now characterize the pairs of wavefunctions at 
energies ±e. Consider running the fermion RG on a given 




FIG. 8. Pictorial representation of the bare distribution 
of the hopping elements: "random walker" steps up by 
ln(Slo/i2n+i) from every odd site 2n + 1 and steps down by 
ln(S7o/t2n+2) from every even site 2n-\-2. Inset shows the ba- 
sic RG transformation when a pair of strongly coupled sites 
2 and 3 is eliminated. Site A survives down to energy scale 
r_B, while site A' is eliminated before this scale is reached. To 
leading order, the A' has component only onto the A, with 
the amplitude given by e~*. 
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system; the concrete disorder realization can viewed as a 
random- walk pattern, while the basic RG step is a trans- 
formation on this pattern that removes the smallest-scale 
structure and keeps the rest of the structure unchanged — 
this is shown in Fig. (note that this is similar to the 
RG approach of Ref. ^ to the problem of random walks 
in random environments). In this picture, the sites that 
remain down to the energy scale are completely spec- 
ified from the initial random-walk pattern as follows: an 
even site A survives down to this energy if and only if 
the "walks" to the right and to the left from the site 
A do not return the same level (of the site A) until af- 
ter they crossed the level F^ higher than the level of the 
site A — see Fig. ^; similar conditions can be formulated 
for odd sites. Moreover, this also gives a corresponding 
criterion for the pairs of sites on which our wavefunc- 
tions at energies ±e are localized (within the zeroth or- 
der RG picture for these wavefunctions). Thus, as far 
as the zeroth order RG picture of these wavefunctions 
is concerned, one is faced with the problem of charac- 
terizing such extremal properties of random walks, and 
the RG itself can be viewed as a powerful tool for doing 
this, even though many questions like the number den- 
sity of remaining sites or the distribution of the effective 
time steps (i.e., bond lengths) and heights (i.e., bond log- 
couplings) can be answered directly in the random-walk 
language by appealing to the properties of a random walk 
near its minima (absorbing walls). 

So far, we have only reformulated our leading (zeroth) 
order result for the wavefunctions at ±e in more visual 
terms. The real usefulness of this picture appears when 
we go beyond the leading order, and keep track of oper- 
ators "measuring" the wavefunction amplitude through- 
out the system. Consider going beyond the leading order 
for the renormalization of the "participation" operators 
rij — When we eliminate sites 2 and 3, the renor- 

malization of the corresponding operators including the 
leading second-order correction to the zeroth order result 
is given as 

(n2)°ff = {ts/t^rui , 
(n3)°^ = (tiA2)'nj. 

It is now easy to see that if at a given stage of the RG 
the site A' is no longer active, it will have component 
only onto the nearest active site A on the same sub- 
lattice, with the magnitude of the "participation" given 
by e~^^, where A is the absolute "height" difference of 
the random walker at "times" corresponding to the two 
sites. The connection to the statistics of the zero en- 
ergy wavefunction is now explicit. Moreover, our picture 
of the wavefunctions at energy e now follows: the dom- 
inant feature of such wavefunctions V'ie are two sites, 
with order one amplitudes of the tp±e, separated by a 
distance of order F^; the wavefunctions will have signifi- 
cant weight only near the two "peaks" and, to this order, 
only on the sites of the same sublattice as the nearby 
peak, i.e., the wavefunction can be thought of as being 



composed of two bumps centered on the two sites. As 
we move away from the two peaks by a distance r ^ F^ , 
the amplitude of the wavefunction is given by the above 
random-walk prescription and is typically decreasing as 
g-c^/r. jTQj. g^jj^j larger distances the amplitudes are given 
by high-order perturbation theory, with the perturbative 
expression for necessarily "stepping" through the 

"active" segments of lengths ~ F^ separating r from the 
two bumps and acquiring a suppression of order e~'^^^ 
from each such segment — thus, for r > F^ we will typi- 
cally have the decay ■0(r) ~ e^'^ ^'^^^ (note that this gives 
us an interpretation, in terms of the typical wavefunction 
at energy e, of the two length scales — one proportional to 
Fe and the other to F^ — that diverge in the limit of low 
e, and are usually identified with the typical and average 
localization length respectively). 
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